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ABSTRACT 


An efficient algorithm has been presented for solving 
the multicomponent, multistage separation process problems 
involving single columns, columns with pump-around or bypass* 
or a system of interlinked columns. In this algorithm I'iiich 
is an extension of the Naphtali~Sandho,lm method, the sparsity 
and the structure of the submatrices in the Jacobian are 
exploited, while performing the matrix multiplications and 
inversions. The operation count was performed for the 
various matrix multiplications and inversions involved in 
the algorithm, and using a variety of tes b-problems it has 
been shown that the sparsity exploitation results in a 
significant reduction in the computational and storage 
requirements. The saving in the computations improves with the 
increase in the number of components. An efficient approach 
has been proposed to solve problems with intermediate tray 
specifications, in which the tridiagonal band structure of 
the Jacobian is preserved, and has been shown to be more 
advantageous in both the computational and storage aspects, 
than the one proposed by Hofeling and Seader; 
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CHAPTER 1 
INTRODUCTION 

The separation processes like distillation, absorption 
and extraction etc, are widely used in the chemical industries 
to fractionate the multicomponent mixtures. Complex systems 
of interlinked columns have been shown to frequently achieve 
more economic and effective separation of multicomponent 
mixtures than the conventional single columns and their 
sequential arrangements. In view of the increasing energy 
and raw material cost, a continuous evaluation of the perform- 
ance is required for economic operation as well as to maintain 
the quality of the product. This requires simulation of 
multicomponent interlinked c olumns. 

The Naphtali-Sandholm method [ 1 ] is widely used for 
the simulation of single columns due to its good convergence 
characteristics. For the simulation of interlinked columns 
a simultaneous rather than a sequential approach is preferred. 
That is» rather than repeatedly solving a sequence of single 
column problems until the solution to the entire system 
is obtained, it is preferred to solve for all the columns 
simultaneously. 

In one version of this approach, the methods like 
capital-theta method of Holland [lo], the equations describing 
each column are considered separately but all these single 
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column problems are converged simultaneously rather than 
repeatedly one at a time. Whereas the more powerful 
approach is to consider the equations describing all the 
separators simultaneously, and solve by the Newton-Raphson 
or a similar method with good convergence characteristics. 
Following this approach several methods have been developed 
by Hofeling and Seader [2], Kubicek [7,12] , Stadtherr [5], 
Browne [13] etc. While these methods differ in their 
formulation of equations and their selection of independent 
variables, all rely on reordering some or all of the linearized 
equations to produce an almost-band or almost-block-tri diagonal 
coefficient matrix. Whereas these methods differ in their 
approach to obtain the desired matrix form and their method of 
solving the reordered linear system^ 

The Naphta li~Sandholm method which was originally 
formulated to solve the single column problems has been 
extended [ 2 ] , to handle the system of interlinked columns, 
and the columns with pump-around or bypass, while retaining 
the technique of total linearization and simultaneous solution 
of all the equations in the system by the Newton-Raphson 
method. 

Though this method offers a good convergence character- 
istics, the large computational and storage requirements are 
the impediments. 

When the Newton-Raphson technique is applied to solve 
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the equations in the Naphtali-Sandholm method, the resulting 
Jacdaian matrix has a tridiagonal~band or a Imost-tri diagonal- 
band structure with a few off-tridiagonal elements. The 
submatrices in the Jac(±)ian are sparse and have a definite 
structure. 

It appears from the literature that the sparsity 
of the submatrices has not been exploited. The objective of 
this work is to present an efficient algorithm, vhich takes 
the advantage of the sparsity and the structure of the 
submatrices in the Jacobian, to solve the separation process 
problems , 

In Chapter 2, first the method of formulation and 
solution for single column problems has been presented and 
then extended to the case of interlinked columns. Chapter 3 
deals with the exploitation of sparsity in the various 
matrix multiplications and inversions. In Chapter 4, an 
efficient approach has been proposed to solve the problems 
with the intermediate tray specifications, and a variety of 
specifications for condensers, reboilers and intermediate 
trays have been presented. Chapter 5 deals with the details of 
implementation of separation process problems on the computer. 
The results of various test-problems and a brief discussion 
have been presented in Chapter 6 and the conclusions in 
Chapter 7^ 
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CHAPTER 2 


MODELLING OF MULTICOMPONENT MULTISTAGE COLIMNS 


In this chapter, first the modelling of a single 
multicomponent multistage separation column is presented, 
which is later extended to a system of interlinked columnsi 
In both the cases, the model proposed by Naphtali and 
Sandholm [1] , is follovjsd, except for a minor reordering 
of variables and discrepancy function in the formulation of 
the , JacdDian-matrix, the reason for which it is done is 
explained later. 


2,1 Generalized Tray Model for a Single Column 

A model for a general tray j, with streams to 
and from the neighbouring trays j-1 and j+l , and with side 
streams is shown in Figure 1 , The discrepancy functions 
formulated using mass balances, equilibrium relations and 
energy balance are as undert 

Component ftess Balances; 


M. . 


= 1 


+ V 


for 1 < i < C and 1 j <. N 

( 2 , 1 ) 
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Equilibrium Relationships: 

Rs K- i i (i “* 

Qj.i = ^3,i - ^ 


Vj+i ''3+1 ,i 


( 2 . 2 ) 


for 1 < i < C and 1 ^ j _< N 


where 


n, = 


^1.1- V 1 .. 1 .. 

''J.i "3,1 - Vi.i 


(2i3) 


Energy Balance: 




0 

®3 ° ii, *’ 3 - 1,1 " 3 - 1,1 " -i, " 3 + 1,1 ■3+1,1 


i A i E 4 4 “* 


i=1 


’j- 


2 h. . 1. . 

i=1 


(1+ sj S 
^ i=1 


H 


c 

2 


:. . V. . + .. h_ f . . + q. 

j,x 1=1 3,1 

(2.4) 

for 1 ^ j jg N 


where, Ij ^ and v. ^ are the molar liquid and vapor flow 
rates leaving the tray j after the liquid and vapor side 
streams s-l^ and S.v. , have been drawn from stage j; 

jj»^ 

Thus there are (2c+1) equations ( discrepancy functions) , 
viz, M4 ^ ,Qj 4 and E., and (2c+l) variables, viz, 1- 4,v. . 

J f J f J- J J f i J f 

and Tj for each tray and hence a total of N(2c+l) of each 
f or the v^ole column. The set of these N( 2c +1 ) discrepancy 
functions and variables may be written in a compact form as 
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follows: 


F( X) = 0 


(2i5) 


where, 


P ~ (F^f ^2* 


( 2 , 6 ) 


and 




(2.7) 


for 1 ^ j ^ N 


X = 

mm ^ mm ••• y 

(x,. Xg, Xj, Xjj) 

(2.8) 

= 

^h.i' h.2 h.c ’ '' 3 , 1 ’ 

'' 3,2 '^3,0 


for 1 < j < N 


(2.9) 


It may be noted here that; 


(a) 


(b) 


The ordering of the discrepancy functions and variables 
is a bit different from the one proposed by Naphtali 
and Sandholm [1] . This reordering results in the 
suitable structure of submat rices for the purpose of 
sparsity exploitation, which is explained later. 


Q. 


A minor change has been iuciu« jui — 

proposed by Naphatali and Sandholm ass 

j+1,i 


'j.i* 

Q. 


which was 

Tl_. 


V_. 






if. 


1 . . - V . . 


+ 




( 2 . 10 ) 
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The discrepancy function (2,2) is obtained by 
dividing the original equation (2,10) by V., This minor 

%J 

change has improved the convergence characteristics. 

Employing the Newton-Raphson method to solve the 
system of equations ( 2,5) simultaneous! y we obtains 


/ dF \ — 1 

I MmCiuMnW* j 

dX m 


( 2,11 ) 



( 2 , 12 ) 


v\iiere, J is the Jacobian matrix and m is the iteration 
number, and, 



The structure of the Jacobian matrix iss 



(2.14) 



8 


where ^ 


3F. _ ' SF, 

■>1 ■ ■ i m rnm it ^ ^ SS I. 


- '^3 = 


(2*15) 


The submatrices , B. and C. have the following elements: 

*J 

3M. . 3M. * 3 M. * 3M. ' 3M. 

*l . ••• 'sn*']'' 3"v^, ••• 

1-1 »1 J-1fC 1-1,1 1-1 


* m M ^ M * M • 

J~1,C 3~1,1 


j-1,c 


®“i c 

^.1 

>1 

j-ul 


3 M . ^ 3 M . 3 M . 

••• rrTT 3 V. ^ 3 v.' 

3-1, c 3-1,1 3~1,C 


9Q. * 3q-. / 

••• 3 1. , ^ vT V . 

3-1. C 3-1.1 


j-1,c 


3M. ^ 




— 3V^ 


• 4 3 V . * . 

3-1. c 3-1.1 


'^3-1.0 


aqv-' 


L%. * . 

3-1,1 


3 3 


3 1, , ^ 3 V . , . 

3-1 .C 3-1.1 


3 V . . 

3-1 .c 




(2.15) 
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CS 


®“1,1 

aM. 


, ^.“l,,1, 

aM ; 7 

ljul 

”3 


••• ^ 

^''3,1 




# 

f 

• 

• 

♦ 

• 

# 

• 

# 

« 

• 

• 

4 



'“i.ro 

3M. 

1«C 

c 

5T^ 



- ^■''3.0 


... ^"i.1 


. ^,'".1,3 




^''3,1 

•*• 3v. 

J»C 

# 

• 

. 1 

# 

# 

« 

• 

• 

41 

• 

# 

• 

• 

' • 

• 

• 


3 Q • 


aq. 

.... — Ll£ 

j j*- 

aq . 

“TT^ 

3 


••• ITT^ 

Jic 

^''3.1 

9E . 

,^1 

a E . 


aE^. 

aE. 


®''3,1 

••• av. " 

J»C 

®"3 


(2w17) 




^"3,,:i: 

*“.1.1 

aM . ■• , 

.1. 1 

® 3+1.1 

4 

gk 

... 

# 

# 

# 

”3+1' 

• 

W 

^Al£ 

« 

... ^ 

• 

• 

♦ 

aM. 

. .1 . c 

• 

• 

aM; - - 

l.C 

^Vi.i 

^,0 

a-M,i 

9 V 

3+1 »C 

® Vi 


.^1*1 


aq. , 

,®^.i,ri : 

3 1-4^ 
J+1 11 

• 

*** 3 1-a.^ 

J+1 fC 

4 

®V1,1 ‘” 
«. 

av . , 

3+1 *c 

^^3+1 

.4 

« 

^1 

♦ . 

3 q’ 

... 

^Vl,c 

• 

aq. 

JlJxE 

9v.. , •■•• 

J-^-t ,1 

• 

a V . . 

3+1 »c 

1 

■f-' 

® c 

Ifi. 

a E^ 

aE. 


dE^- 


••• ‘STTr 

J+1 fC 

^.r” 


■yF± 


( 2 ^1 av 
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In compact form, the structures of A , B and C submatrices 
may be written as 



where, 

is an identity matrix of order cxc 
is a null matrix of order cxc 


( 2 . 19 ) 


( 2 . 20 ) 


( 2 . 21 ) 


6^ is either a row or a column null vector ( whose 


meaning is clear from its position) of order c 
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X is a matrix with non-zero elements, of order cxc 

X is either a row or a column vector with non-zero 

c 

elements, of order c 
X is a non-zero e lament i 

It may be noted that the Jacobian matrix in (2,14) has a 
tridiagonal-band-structure and is highly sparse. This is 
due to the fact that only those derivatives are non-zero 
which are obtained by differentiating the discrepancy functions 
of any tray j, with respect to the variables of the same 
tTay 3 , or the variables of the neighbouring tiays j-1 or 

j+1 » 

Further, the submatrices X and 6 in the Jacobian J 
are also sparse which may be noted in the structures presented 
in (2,19) and (2,21), 

2*2 Interlinked Columns; 

The Naphta li-Sandholm method which was or iginally 
formulated to solve the single column problems, can be 
extended to handle a system of interlinked columns, and the 
columns with pump-around or bypass; 

Generalized Tray Model ; 

A model for a general tray j, in a system of inter- 
linked columns or in the columns with bypass or pump-around 
streams is presented in Figure 2, In addition to the streams 
to and from the neighbouring trays and side-streams, this model 
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includes the interlinked-streams which link the tray j 
to trays which are not its immediate neighbours in the 
system-model. 

Let Rjp Lp, RjqLp, •..be the streams 

v\tiich are leaving the tray p, q, y, z and are fed to 

the stage j, where represents the fraction of stream 

\Atiich is leaving the tray k and fed to tray j, and 



= 0 if there is no stream from tray k to tray j 

= 1 if the whole stream from k is fed to tray j 

-a a value between 0 and 1 if a part of the 

stream leaving the tray k is fed to tray j. 


The sum of all liquid and vapor interlinked streams leaving 
the tray j is represented by (1 - r, )L. . and 
( 1 - r,. )V. , . , respectively. 

The discrepancy functions formulated based on the 
component material balances, equilibrium relationships and 
enthalpy balance ares 


Component Mass Balances s 




r 1 . 




(1+sJl. . ~ (H-S.)v. , 


3 1,1 


3,1 


^3P^P»i ^2,i •••• 


+ fj^i 


for 1 £ i ^ C and 1 < j <, N 


(2422) 
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Equilibrium Relationships: 

Q = 1 _ ^ :: - 

3.1 Lj 3,1 Vj Vy + 


( ^ — T).) R. V ♦ 

V iv v.i 

r,- V. . , + R. V +R. V 
3+1 jy y 3Z Z 


( 1 -n.) R. V ; - 

\ • "1 ^ "t T9^ ^ T 


r,; V.^ . + R.„ V„ + R^^v„ 

3+1 jy y 3z 2 


+ ... (2.23) 


for1 <i<Cand 1 <j<N 


Enthalpy Balances 


"3 = 


y h. . . 1. . . + r,, 
3-1 3-1 V. 


3-1 i=1 ^ ^ ^3+1 


y H.^. . V. 

" l-i-i T -1 


3+1, i 3+1, i 


u u 

-(l+sj E ^-i i 1-i i ~ (1+Si) E H. . V. . 
3 3,1 3,1 3 3,1 3,1 


+ R. V h , 1 . + R. S h_ . 1_ . + 

3P P»^ i=t 


+ R-„ 

3Y 


“y,i ''y.i ^ Hz.i''z.i+ ••” 


+ S hp f . + q. for 1 ^ j ^ N (2^24) 

i=i f'j.i 3,1 3 


The choice of variable remains the same as it was for a 
single column viz, Ij^i* and Tj, in that order. Thus, 
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we have N( 2c+l) equations and N( 2c+1 ) variables which are 
to be solved simultaneously using the Newton— Raphs on method. 

It may be noted here that the discrepancy functions 
of the tray j may also have the variables, of the trays, 

\Atiich are not the immediate neighbours of the tray j, due to 
the presence of int erlinked-streams, bypass or purap-arounds. 
Therefore, the Jacobian matrix may have some off-tridiagonal 
blocks in addition to the tridiagonal’^band. The position 
and structure of these of f-tri diagonal blocks depend upon 
the arrangement of the columns in the system model 
(Hildalgo and Seader [4]), position of the interlinked trays, 
and the direction and type ( liquid or vapor) of the inter- 
linked streams* 

The single stream between two trays (not neighbours) 
may be classified into four basic types; 

(i) liquid stream from an upper tray to a lower tray 

( ii) liquid stream from a lower tray to an upper tray 

(iii) vapor stream from a lower tray to an upper tray 

(iv) vapor stream from an upper tray to a lower tray. 

The reciprocal streams (in which a stream from one 

stage to another st^qe, is matched by a stream of the other phase 

flowing in the opposite direction between the same two 

stages), may be handled as a combination of these four 

types of streams; 
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Consider an arrange nent of interlinked-columns in 
which j and q are the two trays such that j is above q and 
q-j >2, and a stream linking the two. The presence of 
this stream will cause an off~diagonal block to occur in 
the Jacobian matrix, and the position and structure of 
which may be determined considering the four types discussed 
bel ow; 


Type 1 t 


I = 


where 


A liquid stream fi'om the tray j to tray q^ 




Structure; similar to the tri- 
diagonal A blocks 

Position « below the tridiagonal 
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A liquid stream from the 

(Liquid flovdng upward) 



tray q to tray j 



■ * t 

• • ' • r 


'N “N ! 

J 



Structures similar to the 
tri diagonal 
A-b locks 

Positions above the tri- 
diagonal 



A vapor stream Rj„V from tray q to tray j 

J H M 

(Vapor flowing upward) 
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Structure; similar to tridia- 
gonal C-blocks 

Position : above the diagonal 



A vapor stream Vj from tray j to tray q* 
(Vapor flovdng dovjnward) 




Structures similar to tri- 
diagonal c-blocks 

Positions belov/ the 
diagonal 


The reciprocal streams in an arrangement of columns 
cause more than one of f-tri diagonal blocks to occurs in the 
Jacobian, one offdiagonal element per single stream. The 
reciprocal streams may be handled by breaking them up into 
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single streams and the position and structure all off~ 
tridiagonal blocks can be determined by identifying its 
type as discussed earlier. 

Consider a reciprocal stream between two trays » 
a liquid stream flowing from tray q to tray j and , 

a vapor stream flowing from tray j to tray q^ 

The liquid and vapor streams may be identified as type 2 
and type 4 individually, and thus the structure of the 
Jacobian is; 



wh ere. 




Structure: similar to c-blocks 

Position ; below the diagonal 

( corresponds to the vapor 
stream) 


and 
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Structure: similar to 
A-blocks 

Position : above the diagoial 

(Corresponds to the liquid 
stream) 

Thusj the position of an off-diagonal element 
depends only on the direction of the corresponding stream; 

In the Jacobian matrix, it is positioned in the row corres- 
ponding to the tray into v/hich the stream enters, and the 
column corresponds to the tray from which it exits. 

Therefore, if the stream flows upward its position would 
be above the tridiagonal, and below the tridiagonal if the 
stream flows downwards, in a given arrangement of columns; 

The structure of the off-diagonal elements may be 
determined only by knowing v\tiether it is a liquid or a vapor 
stream. An offdiagonal element corresponding to a vapor 
stream has a structure similar to that of a tridiagonal 
C-block, and similar to that of a tridiagonal A-block if it 
is a liquid stream. 

The detailed calculations of the elements of the 
tridiagonal and offdiagonal elements have been presented in 
the Appendix A; 

2,3 Modelling of Condensers and Reboilers 

Although the condensers and reboilers are treated 
like travs, the generalized tray model, discussed in 
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Section 2,1 , cannot directly be applied since in these 
cas^ a few streams are missing and some additional 
specifications are required. In a condenser the liquid 
stream lj__^ , and in a reboiler the vapor stream are 

not present. Moreover, there may be an additional specifi- 
cation for each condenser and reboiler (These maybe specified 
in various ways, 13 of them are included in the Chapter 4)^ 

We assume that there are no interlinking streams 
present in the condensers and reboilers, and the heat-duties 
are specified , (Refer figure 3 for various types of 
condensers and reboilers). 


Partial Condensers; 


Component Ma$s> Balances; 
“j.i = 


1 i C 


(2.25) 


Equilibrium Relationships 


Q. . 


n. K. , V . C1 -’iJ 

— h.i ■ ^ ~v. — ^ ''' 


0 3+1 

1 < i < C 


1+1 

( 2.26) 


Enthalpy Balance; 


E. = r 

1 


c 

S 


»3+i,i Vt.i ^ 


(2.27) 


where, Q- is the condenser-heat-dutyi 
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Total Condensers ; 

In the case of total condensers both the streams 
leaving the condenser are liquid streams, and have the 
same conposition. We use the notation V- for the distillate 

N./ 

(which is liquid). 


Component mass balances: 




(2.28) 


Equilibrium relationships: (Since there is no equilibrium, 

these are replaced by x . . = y. 
equations) ^ 


Q-i i = i - 1* 1 < i < C 


(2i29) 


Enthalpy Balance; 


P . ss r r H V 

J vj+1 j+1,i j+1fi 


I, d.SjVj,! 




4- (l+SjOlj^iJ + 


( 2.30) 


Partial Reboilers: 


Component mass balances: 




(1+S.)v. . , I < i< C 

J J ,X 


(2.31) 
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Equilibrium Relationships: 


K. : 

Q = 1. . 

Lj 3»a. 


V . 

-iii : ■ 

''j ’ 


1 < i < C 


(2;32) 


Enthalpy Balance; 






(2,33) 


where Qj^ is the reboile reheat-duty; 


Total Reboiler ; 

In the case of a total reboiler a part of the 
stream is drawn out as the bottom stream L^* and the remaining 
is vaporized and fed to the column as V., Therefore, the 
bottom stream L. is at the same temperature as that of 
R j_-j > and compositions of Lj and streams are the same. 

Component mass balances: 

“j.i = ^Lj-1 Vi.i ■ ^1^1 = 

(2i34) 

Equilibrium Relationships; (Since there is no equilibrium, 

these are replaced hy x . . = y, . 
equations) J» Jf 


Q . 3= L. V. . - V . 1 . ^ ^ i <. C 

J f -R J J f J J f 


(2i35) 
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Enthalpy Balance: 
c 


i =1 




% 


1+S.) 2 

J i =1 

(2;36) 


There may be various other specifications besides 
the concfenser and reboiler heat-duty specifications, in 
which case the enthalpy balances discussed above will no 
longer be applicable. These are therefore replaced by 
some other suitable equations which take into account 
the given specifications, A variety of specifications and 
alternative discrepancy functions to replace E. are presemted 
in Chapter 4, 


With the replacement of the discrepancy function E . , 

%J 

its partial derivatives will also change. Therefore the 
last rows of the subraatricos (B & C in case of condensers 
and A & B for reboilers) must be modified accordinglyi 


2,4 Method of Solution 

In case of a single column problem with no bypassing 
or pump arounds, the Jacobian has a completel y block- 
tridiagonal structure, and the well knom Thomas-alogorithm 
may be applied to obtain the solution 4 
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Thomas Algorithm; 

Forward Substitution; 

Step 1; ^1 5^ 

Step 2; ** 

For stages 3, from 2 to (N~1) 

Step 3; " ^^3 “ ^3 ^3-1^”^ ^3 

Step 4: Q. * - A. ^3-1^*"^ ^^3 ^ 

For stage N 

Ste p 5 ; ^ % ^N-1 ^ ^ ^-1 ^ 

Ba c kwa r d Subs ti tu ti on ; 

Step 6; 

For stages 3, from CN-1) to 1 
Step 7; AXj ♦ "* 

If in a column bypass or pump-arounds are presen t^ 
or if an arrangement of interlinked column is considered, 
the Jacobian does not have a strictly block-tri diagonal 
structure. Instead, a few of f~tri diagonal blocks occur in 
■,he Jacobian, and consequently, the conventional Thomas 
ilgbrithm cannot directly be applied^ 
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A modification of the Thomas algorithm was developed 
by Hofeling and Seader [ 2], to solve a system in which the . 
Jacobian has a few of f~tri diagonal blocks i 

Hofeling and Seader [2], have taken a specific 
problem to illustrate how the of f-tri diagonal blocks could be 
handled by the modified Thomas algorithm, 

Seader justified that using the same principles of 

the modified Thomas algorithm, any arrangement of the off- 
tri diagonal blocks in a Jacobian could be handled. 

There are a numerous ways, in which the off-tri- 
diagonal blocks may appear in the Jacobian, varying from 
problem to problem. Moreover, for a given system of inter- 
linked columns, there may be a large number of different 
possible arrangement of these off-tridiagonal blocks in 
the Jacobian, depending upon the ordering of the columns or 
column units (refer Hildalgo-Seader [4]) , Although the ways 
to handle these blocks, by the modified Thomas algorithm, 
are the same in principle, it is difficult to develop a 
step by step goieralized algorithm to handle all possible 
arrangements of the of f-tri diagonal blocks in a Jacobian; 

Consider the system of interlinked columns, in which 
there are two absorbers A^ and of 10 plates each, and two 
distillation columns and having 15 and 12 plates 
respectively. These columns are interlinked as shown in 
■-he Figure 6. The Jacobian-matrix, for this system, has six 
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of f-tri diagonal blocks, as shown in Figure 7, 

Vife employ the modified Thomas algorithm to solve 
this system. Various steps of the forward substitution 
and finally of the backward substitution to obtain a 
correction vector are presented here; 


Modified Thomas Algorithm ; 

Step 1 ! row 1 , I 2 q ■** B-| 20* * 1 




step 2; rows j, v\here 2 < j < 9, B. ^ ( B . - A. )"*'*} 

sJ ‘ sJ \J \J ' 

=3' '’3,20 Pj-1,20* «3 *■ ®3^"3-''3°3-1> 


Step 3? row 10, o ^ 0 ^10^9^ 


-1 


10 


^ * ^10,20*^ ~®1 0^1 0 ^9,20* So ** ®lO^So"SoS^ 


Step 4; row 11, P,,,| - B”] ^11,47 ®Tl Si, 47* 


^11,20 ^ ^ 


11 Si 


step 5: row j, where 12 < j ^ 18, B^ '*' ( 


-1 


^ 3 *“ SS’ ^jf20 ^ * ^j,47 “®j^j^3-1 ,47^ 


S* ** ~ 


>tep 6; row 19, ** ^®19 ~ So ^18^ 


-1 
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^19 ■*" ^19^19? 47 ~B.^A 

^ IV 19,47 19^19^18,47* 

^9 - - A , gC !, g ) 

step 7: row 20, set p ^ q . set p „ ► p 

' 2 ^10,20» 


iterate 

on 



iterate 

on 

P2 

- P. 

J 

®20 ** 

(020 - 

^20^19 


^ 20,1 P2 ^ 


-1 


"20 - 5 ; P20.47 *■ - 020 ^ 20 ^ 19 , 47 * 

'^20 ' 020 (020 - A20'5i9 - f^20,^ Pi) 


Step 8: row 21, P„ „ b~| c 


,-1 


21 - “21 ^^21 5 ‘^21 * Bii Fj, ; Pg, 


Step 9! row 22, - Ag^^l )"’ ! 

22 022^22’ 022,47“ '^22 “ 022 (^22 ” ^22^^21) 

Step 10: row 23, - ( Bgg - 


023 * 023^23* 023,47 


®23*23,10 010,20 020,47 


S3 “ 023C 023 - A23S2 “ A23,1o(QiO “ 01 O,2oSo)(> 

Step 11; rows j, where 24 < j < 34 - b . ♦^ ( R - A P 

“ “ * 3 ^3 3 j-1 ^ 


v-1 


“ ®jS= ^j,47 “-Vj^J-1.47* S “Bj (0j-AjQj_P 


tep 1 2; row 35, B 


35 


^®35 ~ ^35^34^ ^ 


35 


^ * 035,47 '- 035 ^ 35034 , 47 * Ss “ ® 35 ( 035 -^ 35534 ) 
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step 13: row 36, , 5 


*- B::3 f 
36 36 ^36 


Step 14s rows j, where 37 < j < 39; B- *(B. - A-P. . 

■J 3 3 3*^* 

Pj ^ BjC^; - 5 ; <5. -B.(F, - 


r 3 J 3-1 


step 1 Ss rov/ 40, B 


40 ■*■ ^®40 “ ■^40^39^ 


reset p ^ •• 

iterate on p ^ ^ j ^ 1 » ^ J’ = 21 

reset ** ^35,47^ 

iterate on ^j,47 “ ^2* 3 “ 34 to j » 23 

p2 **^21^22^2 

P40 ®40^40^ ^40,47'** "’^40^40, 21 ^2’ 

^40 '*"^40^^40 ”■ ^40^39 ~ ^40,21^ 1^ 

Step 16s row 41, 84^ *- (B4,, - A4,,P4 q)''S 

P4I * ®41^41 ’ ^41 ,47 “* ®41 ^^41^40,47'^^41 ,20^20,47^ ^ 

Q4I 341 (^4^ “ ^41^^40 ~ ^41,20^20^ 


Step 175 rows j, where 42 < j < 45, B . ( B . - A .P . . ) 

w J 33 ^ 


-1 




Step 185 row 46, B 


46 


(®46 "" A45 P45) 


-1 


^46 ** ^46 ^^46 *" ^ 46 ^ 45 , 47 ^’ *^46 ® 46 ^ ^ 46 "^ 46 ^ 45 ) 


Step 19 s row 47, ^47 •" ^^4^ — A4 yP 4^) (^^47 “ A4 yQ 4^) 
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Back Substitution; 


Step 1; row 47, ^ 

Step 2; row 46, X,, 4. _ p v 

» 46 ^46 46 47 

Step 3; rows 3, from j = 45 to 3* = 40 

X. ■*-Q.«p,x. -P - Y 
B J J 3+1 ^3,47 ^47 

Step 4: rows 3', from 3 = 39 to 3 = 36 


Step 5; 
Step 6; 

Step 7; 


Step 8s 
Step 9s 






row 35, X^c- 4- Q p Y 
* 35 ^35 35,47 ^47 

rows 3, from 3 = 34 to 3’ = 23 

X.*‘Q.-p.X -p y 
3 3 J 3+1 3,47 "M7 

rows 3, from 3 = 22 and 3 = 21 


6 - ^ ''j - 


^20 ” ^20,47 ^47 


Step 1 0; 
Step 11 ; 


row 20, X20 

rows 3, from 3 =19 to 3 =11 

~ ^J,47 ^47 

ro’" 'O. X^O *Ql 0 - Pl 0,20X20 

rows 3, from 3 = 9 to 3 =1 


***^3 “ ^3 ^ 3+1 “ ^ j ,20 ^20 


It may be noted here that the calculations for eliminating 
the lower of f-tri diagonal blocks is different for the blocks 
^20,1' ^23,10 ■'^40,21 ^41,20* case is discussed 
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here separatelyj 


Row 20 : Block A 2 q^^ ° The structure of Jacobian between rows 1 
and 20 is as under: 



Here, the row of ^ and column of intersect 

on the diagcnal, whereas those of A 2 Q ^ and ^ intersect 
at a point above the diagonal. 

The calculation procedure for this type of situation 
is included in the step 7. Since the C^q block is zero, 

P 2 Q block is also zero. 


Row 23 ; Block 

row 23, is as follows: 


The structure of the Jacobian up to 
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In this case, the row of ^23 10 column of 

20 iritersect at a point below the diagonal, whereas 
those of A23 ^0 Si 40* diagonal. 

The calculation procedure for this kind of situation 
is different from the previous case. Refer step 10, (It may 
be pointed out here that since q and C2Q are zero therefore 
the blocks of rows between 10 and 20 and those of rows between 
20 to 23 do not appear in the calculations) . Please refer 
Figure 7, 

If the blocks q and C2Q were not zeros, then the step 10 
would have been: 


row 23, set P -i •* So* 

iterate on p ^ 3 = 18 to j = IO5 

set p 2 *■ ^19,47* 

iterate on p 2 S 4'7 ~ ^3^2* j = 18 to 3 = II5 

^ 2 ”^^ 10 ^ 2 ’ 
set P 3 ■*' ^ 19 , 20 

iterate on P3 V* ^3,20 ~ ^3 ^2* ^^^3 = 'I 3 = 11 ? 

^ 3 ** ^ 10 ^ 3 * 

set a ^ ^ Q22 

iterate on a ^ - Pj , from 3 = 21 to 3 = 20 

set a 2 ♦* ^22,47 

iterate on a2 **^3,47 ~ ^3 °^2* 3 = 21 to 3 = 20 

set a3^ P22 
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iterate on ^ ^ j = 21 to j = 20 

B23 - (B23 - A23P22 ~ A 23 ^^Qp 3 a 3 )“'' 

^23 *“ ® 23 ^ 23 » 

^23,47 ^ ®23[ ^23,10 ^ ^2 “^3*2^ ~ ^23^22,47 ^ 

^23 ** ®23 [ ^23 “ ^ 23^22 ~ '^23,1 0 ^ ^1 

Thus, the calculation procedure for this type of 
situation, becomes some what more complicated than the one 
discussed earlier. 


Row 40 ; Block A^q 2 j ° structure of the Jacobian upto 

row 40 is? 



MO, 21 


Here, the row of A^q 21 column of 

intersect at a point above the diagonal and there are no 
other of f-tri diagonal blocks betwen the rows 21 and 39 
after the forward substitution. Also the block C 21 is 
non-zero. 
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The calculation procedure for this situation is 
described in the step 15^ 

Row 41 ; Block 2 q* 

Since this situation is similar to that of the row 40, 
the calculation procedure followed is the same. However, 
it gets simplified to step 16 since the block C 2 Q is zero, 
and therefore the blocks between rows 20 and 40 do not enter 
in the computations. 

Thus, the numerous possible arrangements of off-tri- 
diagcnal blocks v\tiich require different computational efforts, 
and the presence of some null blocks which may change the 
computations substantially, make it extremely difficult to 
develop a ’universal* algorithm which can handle all possible 
arrangements. 
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CHAPTER 3 

SPARSITY EXPLOITATION 

The conventional and modified Thomas algorithms 
exploit the sparsity of the Jacobian matrix which is highly 
sparse, and has a completely or almost~tridiagonal-band 
structure , 

Moreover, the submatrices within the Jacobian are also 
sparse (Particularly A and C), and have a definite structure 
(refer 2,19, 2,20 and 2,21), The presence of zero and identity 
blocks within the submatrices and their definite structures, 
allov/ a further exploitation of the sparsity, \Ahich conse- 
quently results in a substantial reduction of both ccmputational 
and storage requirements. 

In this chapter, the exploitation of sparsity within 
the submatrices along with saving in storage and operation- 
count, is presented. 

As is customary, only the operations involving multi- 
plications and divisions are counted, neglecting those 
involving addition and subtraction since the time consumed 
in these is comparatively negligible. 

To take the advantage of the sparsity and structure of 
the submatrices, ccmputati on and storage involving zero and 
unity blocks is avoided, viiich results in a substantial 
reduction in the CPU-time and memory-requirement on a computer. 
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Listed below are the various operations involving 
multiplication and inversion of submatrices, while exploiting 
their sparsity and structures, and comparison, of the operation 
counts. 


It may be pointed out that the P-blocks have two 
different structures viz, Py and Pj^, which are as under; 




fo. X_ 


K 5 ' 

o 

O 

o 


c c c 

0. X_ X 

* p" = 


o 

o 

o 


c c c 

5 x^ X 


X S X 

L c 


c c 


I ~1 / Vs~1 

Also, B = B or (.B-A P ) , where both have the same 


structure 


X, 


X. 


X 

X 


/V 


Operation 1 ; 


B C 


1 X 

ibx 

lie 

\ 




0^ r 5^' 

o 

p 

) 


c c c 


c c c 

^ K ^ 

4 * 

K 


5^ 5^ 

^c c c 


c c c 


c c c 

o X X 


^r. K ^ 



( 

o 

o 


C O 


c c 

L J 

1 





i , k+c 


^\,k 


+ 


2c+1 

S 

m= c+1 




C , 

m m, k+c 


1 ^ i ^ 2c'i-l 
1 < k< c 

waari' tiaa0 
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V 




^i , 2c+1 


^2c+1 , 

,2c+1 

s.o.c. 

= 8C^ + 12C^ 

+ 6C + 

1 

I.O.C.+ 

= 2C^ + 5C^ 

+ 4C + 

1 

Saving 

= 6C^ + 7C^ + 

2C 



1 < i < 20+1 


Operation 2i 



B^A 


K 0^ 

c c c 


c c c 


H 

Mi 

o 

Oil 

o 

oi 

o 

1 

K 


X f X 


on 

on 

OI 

c c c 


c c c 


c c c 



i 5 ^^ X X 


X 0 X 

c c 


c c 


c c 


^i,k+c ** 
^i,2c+1 


^ ®i,k ^i,2c+1 ^2c+1,k 
A 

i,2c+1 ^2c+1,2c+1 


1 <i<2c+1 
1<i^c+l 


5.0. C. = 8C^ + 12C^ + 6C + 1 

1.0. C. = 4C^- + 4C + 1 

Saving = 8C^ + 8C^ + 2C 


♦ Standard Operation Count 
+ Improved Operation Count, 
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Operation 3: R *• A 



1 < k < c+1 

V ^ V 

^2c+l ,k+c *■ ^2c+1 ,2c-i-1 ^2c+1 ,k+c ^2c+1 ,m^m,k+c 

1 < k < c+1 

HaM tstit 

5.0. C. = 8C^ + 12C^ + 6C + 1 

1.0. C. = 2C^ + 2C + 1 

Saving = 8C^ + 1 OC^ + 4C 
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^2c+1,k ^2c+1 ,2c+1^2c+1 ,k ^2c+1 ,m ^m,k 


1 < k< c 


L ^ 

^2c+1,2c+1 "^20+1,20+1 ^2c+1,2c+1 "^20+1,111^^,20+1 

5.0. C. = 8C^ + 12C^ + 6C -!- 1 

1.0. C. = 2C^ + 2C + 1 

Saving = 8C^ + 1 OC^ + 4C 


Operation 5: R ♦- C 





5.0. C. = 8C^ + 12C^ + 6C' + 1 

1.0. C. = C^ + 3C^ + 3C 4- 1 
Saving: = 7C^ + 9C^ + 3C 
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1 , 1 c i+c,k 

^i,2c+1 ** ^ ^i-^c,2c+1 


1<i< c, 1 <k<c 


1 ^ i ^ c 


L 

R.^ 1 2 C.^ P . 1 < i < c, 1 < k < c 

i+c,k i+c,m m,k - - * - - 


\-!-c,2c+1 "^j^c+1 ^i+c,m ^m,2c+1 ^ 


2c+1 , k 


2c+1 T 

? , / P ,, 1 < k < c 

^ ,t 2c+1 ,m m, k — — 

m=c-i-1 * 


2c+1 L 

^2c+1 , 2c-i-1 ^2 c-!- 1 ,m ^m, 2c+1 


5.0. C. = 8C^ - 1 - 1 2C^ + 6C + 1 

1.0. C. = -I- 3C^ + 3C + 1 

Saving = 7C^ + 9C^ + 3C 
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Operation 7; P B-'-R, where R = -AP^ (refer OP 3) 



VI ^ T 

“P ^ R"*" R 1 yi R*^ R 

i,k+c i,2c+1 2c+1,k+c °i,in m, k+c 

1 < i < 2c+1 ,1 <k<c4-l 

«aEi^ MhM * tawa< taaM 

5.0. C. = 8C^ + 120^ + 6C + 1 

1.0. C. = 2C^ + 5C^ + 4C + 1 

Saving = 6C^ + 7C^ + 2C 


Operation 8: 


P^ - B^R, where R = 


A P^ (refer OP 4) 



P^ 

^ i, 2c+1 


‘®i,2d+1 ^2c+1,2c+1 


\,ra ^m,2c+1 


■k i<2c+1 


5.0. C, = 8C^ + 12C^ -i- 6C+1 

3 2 

1.0. C. = 2C + 5C + 4C+1 

Saving = 6C^ + 7C^ + 2C 
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V T 

Operation 9: P i»- B R , where R = C 



5.0. C. = 8C^ + 12C^ -s- 6C + 1 

1.0. C. = 4C^ + 8C^ + 5C + 1 

8 o 

Saving = 40'" + 4C + C 


Operation 10; - B^R where R = C P 



L I 

^i,2c+1 ^ ®i,ni V 2 C- 1 -I 


5.0. C. = 80^ + 12C^ + 60 + 1 

1.0. C. = 4C^ + 8C^ + 5C + 1 

Saving = 4C^ + 4C^ + C 


P^ (refer OP 5) 



, 1<k<c+1 


^ (refer OP 6) 



2c+1 , 1 ^5 k 5. c 

^ i < 2c+1 
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Operation 11 : R 4- pL 



R. 2c+1 V _L 

i,2c+1 *- Z i,mV2c+1 

m=c+l 

S. O.C. = 8C^ + 12C^ + 6C -I- 1 
I.O.C, = 20^ + 5C^ + 4C -1-1 
Saving = 6C^ -t- 7C^ -i- 2C 

Operation 12; R 



R pi-j pJL L L» 

i-,2c-f1 **1,20-1-1 2c-H,2c-M + 2 ^i,m ^m,2c-H ^^-5^0+1 

m=1 

5.0. c. = 80^ -f 12C^ -f 6C -I- 1 

1.0. C. = 2C^ -f 5C^' 4C + 1 
Saving = 6C^ -i- 7C^ -i- 2C 
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Operation 13; R 



1<k<c+1 

S.O.C. = 8C^ + 12C^ + 6C + 1 
^•0.0, = 2C3 -h 5C^ + 4C + 1 
Saving = 6C^ + 7C^ + 2C 


Operation 14; R •*- 



m = c+1 


5.0. C. = 80^ + 12C^ + 6C+1 

1.0. C. = 2C^ + 5C^ + 4C + 1 
Saving = 60^ + 7C^ + 2C 
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Operation 15s S aQ 



5.0. C, = 4C^+ 4C + 1 

1.0. C. = 2C -I- 1 

Saving = 4C^ + 2C 


Operation 162 S <*• C Q 



®i-^c 


2c 

2 

m=c+1 


i+c,m 


1 ^ c 


^ 2c+1 ,m 

m=c+1 

5.0. C. = 4C^ -I- 4C;- 1 

1.0, G, = C -f 20 + 1 

Saving = 3C^ + 2C 

■rnii. 
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Operation 17: S •* Q 



S,0*C, — 4C + 4C + 1 

I,0,C, = 2C 4- 3C + 1 
o 

Saving = 2C +0 


Operation 18: S'** P^Q 



5.0. C. = 4C^ 4- 4c -i- 1 

1.0. C. = 2C^ 4- 3C + 1 

Saving = 2C^ 4- C 
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Operation 19; S - B^F (Standard Operation) 



5.0. C, = 4C^ + 40^ + 1 

1.0. C. = 40^ + 40^ + 1 

Saving = 0 


Operation 20; Inversion of B or (B ~ AP) 

The B and (B~ AP) matrices have the same structure, 
which is as under; 



Inversion by partitioning does not result in saving of 
the computational effort. It may be noted here that this 
matrix contains a cxc identity matrix, I„, and hence some 
saving in the operations may be achieved as demonstrated below. 

We partition the matrix as shown above ( one parti on 


only) 
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Hence, 



or (B AP)-'' 



w^ers , 




I 

I 

( 

{ 

T 

( 

! 

\ 



or der of 
order of 
order of ^2^ 
order of B22 


= cxc 
= cxCofl) 

= ( C+1 ) X c 

= ( c+ 1 ) X (c+ 1 ) 


and| 


^22 = 

^522-521 (E„)-l B, 2)-’ 

(1) 

Bi, = 

(B,,)-’ MB,,)-’ 5,262262, (Bj,)-’ 

(2) 

Bi2 = 

- '-Sn)"’ h2 S 22 

(3) 

^21 = 

-®22 ®21 ^ h 1 ) ' 

(4) 


Since (b^^) = ~ I^.,( only one division is 

needed to compute it), we need to perform only one inversion 
of the (c+1) X (c+1) in equation (1), to compute the 


5.0. C. = 8C^ + 12C^ + 6C + 1 

1.0. C. = 5C^ + 12C^ + 8C + 2 

Saving = 3C^ - (2cH-1 ) 

JS 3C^ 

In the case of a single column problem with no pump- 
around or bypass streams, only the operations 1,3,15,17,19 

and 20 are performed in the conventional Thomas algorithmi 
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Moreover, the fractions r^^ and r^^^ are unity for all the 
trays, since there are no interlinked~streams present. 
Therefore, in the operations 1,3 and 15 there is no need to 
multiply these fractions and a further reduction in the 
operation count is obtained. The improved operation count 
for these operations is presented below: 


Operation 1 ; 
Operation 3 : 
Operation 15 ; 


Improved operation count 

20^ + 3C^ + 3C + 1 
+ 20 + 1 
C + 1 


Saving in Storage Requirement ; 

In the conventional or modified Thomas algorithm, 
computations are performed stage by stage. Since the 
computations for any row j,. require only tri diagonal and 

off-diagonal P-blocks of the previous j-1 rows, the A, B and C 

♦ 

submatrices need not be stored, nor the identity blocks on 
the diagaial after the forward substitution. Thus, we 
need to store only the P-blocks, Since in both types of 
P-blocks, there are ( 2c+1 ) x c zero elements which need not be 
stored. 


r°c 
i % 

t 

‘ 0 
i C 






X. 


■_ 1 _ I 

: 5c I 

: 8c ' 

1 I 


5? 


C 1. °c ^ 




X 
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Thus, out of C 2c+1 ) X (2c+1) elements of any P-block 
only (2c+1) x ( c+1 ) are stored, and thereby almost a 50°^ 
further saving in storage results (Refer Chapter 6), 
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CHAPTER 4 
SPECIFICATIONS 

A multicomponent, multistage separation process 
problem must have the follovdng specifications; 

( i) Number of trays 
(ii) Number of components 
( iii ) Column pressure 

(iv) Complete specification of the feed (the feed rate, 

composition, thermal conditions, and the location of 
the feed tray) , 

In case of distillation, the location and type of every 
condenser and re boiler must also be specified. 

In addition to these, some more specifications are 
needed to completely define the problem. These specifications 
may be made for the condensers, reboilers or for any other 
intermediate stages* It is important that the total number 
of these specifications must be equal to the total number of 
condensers and reboilers present in the system, otherwise 
the problem would become under or overspecified. 

Specifications for Condensers and Reboilers ; 

In the Naphta li-Sandholm model > the specifications 
for condensers and reboilers are their heat duties, which are 
incorporated in the enthalpy balance discrepancy function, E^, 
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For specifications other than the heat duties, the enthalpy 
balance in equations (2.27), (2.30), (2,33) or (2,36) shall 
no longer be applicable, and should therefore be replaced 
by the suitable discrepancy functions which takes into 
account the given specifications. 

Some of these specifications and their corresponding 
discrepancy functions are presented in this chapter. 

Specifications for the Intermediate Trays 

When specifications are made at the intermediate 
stages other than condensers and reboilers, the heat- 
duties are no longer variable, and get fixed accordingly. 
Consequently, the enthalpy balance discrepancy functions, 

Ej, for the condensers and reboilers are no longer applicable, 

Hofeling and Seader [2], have proposed that the inappli- 
cable discrepancy function E.. , for a condenser or reboiler 

J 

tray, maybe replaced by a suitable discrepancy function, 

\Ajhich incorporates the specification at the given intennediate 
stage. This replacement however gaierates' an of f-tri diagonal 
block in the Jacobian since the function for a condenser 
or reboiler contains the variable(s) of an intermediate stage. 
The modified Thomas algorithm may be employed to handle 
the of f-tri diagonal blocks, 

.^E.NITRAL 

.. ^ Kanpur, 
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Generation of the of f-tri diagonal blocks results 
in additional computational and storage requirements, and 
should therefore be avoided, if possible. It is therefore 
proposed here, that the discrepancy function E- for condenser 
or reboiler may be modified and made applicable (as in 
Spc. 1 and Spc, 2 below), and the E. for the intermediate 
stage should be replaced by a suitable discrepancy function 
to incorporate the given specification at the same tray. 

Listed below are som.e specifications and corresponding 
discrepancy functions which may be used to replace the 
inapplicable ones: 

Spc, 1s If a specification is made at an intermediate 

stage instead of condenser, the condenser heat duty 
Q in equations 2,27 and 2,30 cannot be specified and 
should therefore be computed. Calculations of Q^, 
and the modified are as under; 

Partial Condenser; In a partial condenser, (which is 
represented by tray j), out of the r^^^^ , 
moles of the vapor stream, (1 + Sj)Lj moles are 
condensed (Figure 3). Therefore, maybe computed 
by enthalpy considerations s 
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Q, 


and 


=3 = 




il, Vl.i ^"vj+1 

c 

~ (1 + s.) 2 H. . V. . 

1=1 J 5 ^ J f ^ 


Total condenser; In a total condenser, tray j, the wAiole 
^vj+i ^j+1 stream is condensed. 


Q. 


=3 


■^vj+1 ^^3+1.i " '^3.i^''3+1.i 


, i ^ ^3-1 "3+1 . 1 - ‘ '^=3 ^ ^3 , 1 ' ^3 ’"3 .1 ’ 


Spc, 2; If a specification is made at an intermediate 
stage instead of a reboiler, the reboiler heat duty 
Qj^ in equation (2i33) and (2.36), cannot be specified 
and should there fere be computed. 


Partial Reboiler; In a partial reboiler (which is 
represented by tray j), out of the moles 

of the liquid stream, (1 + Sj)Vj moles are vaporized 
(Figure 3). may be computed by enthalpy consider- 
ations ; 


Qr = Cl 


+ Sj) 


C 

2 

i=1 


( H . • — h . . . )v . . 
^ j,i 3~1»i 
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"r 


c 

vj-i 3-1,1 3^ 3,1^ 

c 

- (1 -i- s^) 2 h. 1. 

izr-t J J f * 

Total Reboiler; In this case, out of r . ^ L . ’ 

' vj-1 3-1 

moles, (l + Sj)lj moles are vdthdrawi and remaining 
(1 + Sj)V^. are fed to the re boiler and are completely 
vaporized (Figure 3), 


Q 


■R 




^0 h-i.i ■ 

Spc, 3; Specification of the condenser-heat-duty, Q , 

C 

As discussed earlier, (Please refer equations 
C 2,27) and (2.30) for E^). 

Spc, 4; Specification of the reboiler-heat-duty, Qj^, 

As discussed earlier (Please refer equations (2^33) 
and ( 2,36) for E .) 

sJ 

Spc, 5; Specification of the Reflux Ratio, R = L/D for a 
condenser. ( or R = (L/V) may be specified for an 
intermediate stage) 

Ej = L . - R V. 

Spc, 6s Specification of the Reboil Ratio, r = V/B for a 
reboiler. ( or r = (V/L) may be specified for an 
intermediate stage) ^ 
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Spc , 7; SpGcixica'ti on of Teraporature at any stags j 

E. = T . - T 
J 3 spec. 

Spc, 8- Specification of total vapor flow rate at any 
stage j. 

^^spec ~ ^ “ Distillate rate, for a condenser) 

E . = V- - V 
3 3 spec 

Spc, 9; Specification of total liquid flow rate at any stage j, 

= B = Bottoms rate, for a reboiler) 


' spec 
E^ = L, - L 


spec* 


Spc, 10; Specification of molar vapor flow rate of component 

i at a stage *j. Vj,i)spec. = ^i = vapor key 
component flow rate, for a condenser tray), 

^3 = ''3,1 - ‘''3,i5spec- 

Spc.11; Specification of molar liquid flow rate of component 


i at a stage j. 


liquid key 


component flow rate, fora reboiler tray), 

^3 ~ ^3»i ~ ^^3»i^spec* 

Spc,12; Specification of vapor mole fraction of component i 

at stage 3. ((yj,i)spec. = ’'d. = 
component composition for a condenser tray), 

p.5=v..~V.(v--) 

^3 3»i 3^^3»i^spec. 

Spc. 13; Specification of liquid mole fraction of component i 
at stage j. "" ^b^ "" liquid key component 

composition for a reboiler tray) 

^i,i " ^3 ^^3,i^sp€C* 
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Consider a single distillation column vdth 15 
trays (including the condenser as tray 1 and reboiler as tray 
15), and a 6-component mixture. Instead of condenser and 
reboiler heat-duties, the temperature of tray 9 and the flow 
rate of vapor leaving the tray 10 are specified,. 

As proposed by Hofeling and Seader in [2] , the 

condenser enthalpy balance discrepancy function , is replaced 

by E. = - (Tf^) . ,, and E. p. of reboiler by 

^1 9 9^ specified* 15 * 

^15 “ ^10 “ ^^1 0^ specified* resulting Jacobian structure 

is shown in Figure 8, which may be solved using the modified 
Thomas algorithm, 

Altematively,we can modify the discrepancy functions 
and as discussed in Spc. 1 and Spc, 2, and replace those 
of rows 9 and 10 by E, = Tg - and = 

Vi 0~^Vi q) specified’ ^Qspectively, The resulting Jacobian 
matrix, in this case, has a completely tri diagonal structure 
with no of f-tri diagonal elements. The conventional Thomas 
algorithm is employed here, and a substantial reduction in 
the computational efforts (3,877 times), and in memory 
requirements ( 2,786 times) is thereby realized. 
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CHAPTER 5 

IMPLEME.NTATION OF THE ALGORITHM 

A computier program for solving the multicomponent 
separation process problems .has been written in FCRTRAH 10 
and implemented on DEC 1090 system based on the algorithm 
described earlier. The program is kept in a file SEP, FOR, 

The file SEP. FOR consists of a main program INTLNK 
and various subroutines and the tasks performed by each of them 
are listed belows 

Main Program s 

INTLNK; The program may be considered to have two parts 
with different tasks. The first part reads the 
data, and computes the discrepancy function - vector. 
This part has been written for a general system of 
interlinked-columns and can handle any arbitrary 
arrangement of columns. 

The second part includes the algorithm for solving a 
specific problem. As has been pointed out in 
Chapter 2, it is difficult to develop a generalized 
Hofeling-Seader algorithm which can handle all 
possible arrangements of the off-tridiagonal blocks in 
a Jacobian, Therefore this part must be rewritten 
for a given arrangerrent of the off-tridiagcnal blocks. 



58 


Subroutines 

ENL ; Computes the liquid— molar-enthalpy for all 

components at a given stage (using a polynomial), 
ENV ; Computes the vapor-molar-enthalpy for all 

components at a given stage (using a polynomial), 
FINDK ; Computes the equilibrium values for all 

components for a given stage (using Antoine’s 
constants (ideal-case)); 

DKBYDX ; Computes the derivatives 

a X 

3 k 

DKBYDY ; Computes the derivatives ^ 

a y 

3 If 

DKBYDT ; Computes the derivatives ^ 

0 H 

DHVDT ; Computes the derivatives ^ 

DHLDT ; Computes the derivatives 

CONBC I Computes the elements of the B & C submatrices 

for a par'tial or a total condenser, 

REBAB ; Computes the elements of the A & B submatrices 
for a partial or a total reboiler, 

TRIA ; Computes the elements of the A submatrix for a 
given tray, 

TRIB s Computes the elements of the B submatrix for a 
given tray, 

TRIG ; Computes the elements of the C submatrix for a 


given tray, 
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OFFA s Computes the elements of the off~tridiagonal 
block having the same structure as that of an 
A-submatrix, This block is generated because 
of the liquid in ter linked- stream, 

OFFC ; Computes the elements of the of f-tri diagonal 
blocks having the same structure as that of 
C-matrix, This block is generated because of 
the vapor -interlinked stream. 

PROS i Replaces the last row of the A, B and C sub- 
matrices to incorporate the given specification 
for a stage j, where j may represent a condenser, 
reboiler cr an intermediate stage. 

The follomng subroutines perform the various operations', 
which have been presented in the Chapter), 


BCMUL 

o 

o 

Operation 

1 

BAMUL 

o 

o 

Ope rati on 

2 

APmUL 

■o 

o 

Operation 

3 

APLMUL 

o 

o 

Operation 

4 

CPVMUL 

o 

o 

Operation 

5 

CPLMUL 

o 

o 

Operation 

6 

MBAPV 

0 

Operati on 

7 

MBAPL 

o 

o 

Operation 

8 

BCPV 

D 

0 

Operation 

9 
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BCPL 

O 

O 

Operation 

10 

PVPL 

o 

o 

Operation 

1 1 

PLPL 

o 

a 

Operation 

12 

PLPV 

6 

* 

Operation 

13 

PVPV 


Operation 

14 

AQMUL 

o 

Operation 

15 

CQMUL 

o 

o 

Operation 

16 

QMPVX 

o 

o 

Operation 

17 

QMPLX 

o 

Operation 

18 

BFMUL 

o 

Operation 

19 

INVPRT 

c 

0 

Operation 

20 (This subroutine computes 



the inverse of a B or (B-AP) matrix 



of order ( 2c+1 ) x ( 2c+1 ) ) 

FiATINV 

o 

0 

Computes 

the inverse of ( c+l )x( c+1 ) 



matrix using Gauss Jordan method 



with the 

maximum pivot strategy 



(called by INVPRT). 

After the 

input- 

•data has been read, a check is 


conducted to find whether the problem is under or over 
specified, and if found so the execution is immediately 
terminated with an error message. Several checks are made at 
the various stages of execution, and if any inconsistency 
is detected, a suitable error message is printed out and the 
execution is stopped. 
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For a single column prdDlem, the variable NOFF, 
which represents the number of of f~tri diagonal blocks in 
the Jacobian, should be set to zero, and the program 
uses the conventional Thomas algorithm to solve the 
problem. 

The thirteen different specifications for the 
condensers, reboilers and intermediate trays, which are 
presented in the Chapter 4, have been implemented 
in the program. For a given specification, a suitable 
code (l through 13) is assigned to the variable LABEL and 
the specifications are handled following the approach 
discussed earlier. The subroutine PROB modifies the corres~ 
ponding derivatives in the last row of the submatrices A, 

B, and C to incorporate the given specificationw 

In a system of interlinked columns a few vapor 
or liquid streams may be missing. This situation may 
arise in the cases viflsere splitting and rearrangement of 
columns is done to obtain the *'best~ordsring*' [4] . The . 
locations of the trays for which a leaving-vapor stream 
is missing are stored in the JNOSV array, and those 
for the missing liquid streams in the JWOSL array. The missing 
streams are omitted in the calculations. 

The fractions r^ and r^^ (refer Chapter 2), are unity 
for the trays from vhich no interlinked stream is leaving. 
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an(® for all single column problems* Therefore these 
need not be included in the various multiplication 
operations such as operation 1,2 and 3 etc. The subroutines 
take care of this fact and a further saving in the operation 
count may be realized. 
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CHAPTER 6 

RESULTS AND DISCUSSIONS 

The exploitation of sparsity of the submatrices, 
while performing the matrix multiplications and inversions 
in the algorithms presented in Chapter 2, has resulted in 
a significant reduction in computational and storage require- 
ments, The reduction in both, the operation count and the 
storage, for a variety of problems is presented in this 
chapter. 

Problem 1 ; Consider a single absorption column having 
8 stages and using a mixture of 14-components (For the 
complete specifications please refer [93), 

The problem was solved by applying the conventional 
Thomas algorithm, and a comparison of the operation counts 
is made in the Table 1 , 

A substantial reduction in the computation by a 
ratio of 3,02 was obtained by the sparsity exploitation, as 
only 183774 operations were performed instead of 555060, 

Since only 3045 (= 7x29x14) elements need to be stored 
instead of 5887 (= 7x29x29), a saving of about 4^ in 
the storage was realized. 

This problem when implemented on DEC-10 took 

-4 

4 iterations to converge to a tolerance limit of 9x10 
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(sum of squares of all the discrepancy functions). The 
actual CPU time spent in the matrix multiplications and 
inversions v^ere 32.59 and 11.61 seconds for the standard 
and improved operations respectively. Thus the ratio of 
the CPU time is about 2.81 as compared to the theoretically 
obtained ratio of 3.02, 

The total CPU time spent in the execution of the 
programs for the standard and improved operations are 38.43 
and 17.43 seconds respectively, and the balance was spent 
in evaluating the discrepancy functions, elements of the 
submatrices, and in executing the various input and output 
statements. 

VJith the original discrepancy functions, proposed by 
Naphtali and Sandholm (refer equation 2,10), this program 
took 7 iterations to converge, whereas only 4 iterations were 
needed vdth the modified equation 2,2. Thus, this minor change 
has improved the convergence characteristics. 

It was realized that the large roundoff errors in 
the inversion of the submatrices may impair the convergence 
characteristics. Three different subroutines used to study 
the convergence characteristics, and the results of which 


are as under; 
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Matrix inversion subroutine Number of iterations 

r eqiHFedT 

1 . Gauss Jordan with the maximum 4 

pivot-strategy 

2* Gauss elimination 5 

G, FOlAAF/NAG (available on DEC-10) 12 

A check on the roundoff errors revealed that the 
FOl AFF subroutine resulted in substantially large roundoff 
errors, in the matrix inversions. 

It was also noted that the inversion by partitioning 
method resulted in less roundoff errors as compared to the 
standard matrix inversion, particularly in the matrices where 
the difference in the magnitude of the elements was vast. 

Problem 2 ; Consider a single absorption column having 
20 stages with a mixture of 4 components, (For the complete 
specifications please refer [1J), 

This problem was solved using the conventional 
Thomas algorithm and a comparison of the operation counts 
is made in the Table 2; 

A reduction in the computations by a ratio of 2,676 
is obtained by sparsity exploitation, as only 17556 operations 
were performed instead of 46980, 

This problem when implemented on DEC~10, took 

-8 

4 iterations to converge to a tolerance limit of 5x10 
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TABLE 2 


Comparison of Operation Counts for Problem 2 



Number of 

Number of 

components = 4 

stages = 20 


Operations 

Number of 

Standard 

Improve d 


times 

operation 

operation 


performed 

count 

count 

Op 1 

19 

x729 = 13851 

x189 = 3591 

Op 3 

19 

x729 = 13851 

X 25 = 475 

Op 15 

19 

X 81 = 1539 

X 5 = 95 

Op 17 

19 

X 81 = 1539 

X 45 = 855 

Op 19 

20 

X 81 = 1620 

X 81 = 1620 

Op 20 

20 

X 729 = 14580 

x546 = 10920 

Total 


46980 

17556 


5.0. C . 

1.0. C. 


2,676 
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(sum of squares of all discrepancy functions). The total 
CPU time spent in the matrix multiplications and inversions 
were 2,82 and 1,22 sec, for the standard and improved 
operations respectively. Thus the ratio of the CPU time 
is about 2,31 as against the theoretically obtained ratio 
of 2.676, 

The total CPU time for the standard and improved 
programs was 7,02 and 4,34 seconds respectively, and the 
balance was spent in the evaluation of the discrepancy 
functions, the elements of the submatrices, and in executing 
the input and output statements. 

There are a total of 19 P-blocks in the matrix 
after the forward substitution, each of order 9x9, Since 
only 19x9x5 elements of the P-blocks need to be stored 
instead of 19x9x9 4A% saving in the storage is obtained. 

Problem 3 ; Consider a single column with 15 stages and a 
6~component mixture. In this case, as discussed in the 
Chapter 4, instead of specifying the heat duties of the 
condenser and reboiler, specifications at the 9th and 10th 
stages were made. The temperature of the stage 9, and the 
flow rate of the vapor leaving the stage 10 are specified. 

In the Chapter 4, two alternate approaches to solve this 
problem have been discussed, and the operation counts for 
both are compared in the Table 34 
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For the first approach which was proposed by 
Hofeling and Seader [2], the standard operation' count is 
computed, for obtaining the solution by the modified 
Thomas algorithm using the Jacobian structumre shown in 
the Figure ^ 

For the latter approach, the improved operation 
count with the sparsity exploitation was canputed for 
obtaining the solution by the conventional Thomas algorithm. 

The Jacobian matrix in this case has a tridiagonal band 
structure. 

The standard and improved operation counts are 
145509 and 35729 respectively, v\hich results in a computation 
reduction ratio of 3,877, 

In the Hofeling-Seader approach, there are 21 
P~blocks (14 tri diagonal and 7 off-tridiagonal), each of 
order 13x13, Whereas in the second approach there are no 
off-diagonal blocks and hence a total of only 14 P-blocks in 
Jacobian, 

Therefore follovdng the latter approach, only 
1274 (= 14x13x7) elements need to be stored, whereas the 
standard storage requirement for the two approaches is 
3549 (= 21x13x13) and 2366 (= 14x13x13) respectively. 

The proposed approach is therefore more efficient, 
then the Hofeling-Seader approach since it requires substantially 
less computations and storage. 



Bi Cl 
A 2 Bg 
A 3 


C 1.9 
Ca 

B 3 Ca 
A 4 B 4 C 4 

As Bs Cs 
Ag Bg Cg 

Ay B 7 C 7 

Ag Bg Cg 

Ag Bg Cg 

Aid Bio Cio 

An Bn Cji 

Ai2 Bi2 Ci2 

Ai3 Bi3 Ci3 

Ai4 Bi4 Ci4 

A 15.10 Ai 5 Bis 


FIG 4 THE STRUCTURE OF THE JACOBIAN MATRIX 
(PROBLEM 3) 
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Problem 4 ; Consider the system of interlinked columns 
shown in the Figure g", in which a 4-component mixture 
is being separated using two absorbed and A 2 having 
10 plates each, and two distillation columns having 15 and 
12 plates respectively. The arrangement of columns and 
ordering of the trays is also shown in the same figure, 
and the structure of the Jacobian matrix is presented in 
the Figure (For the complete specifications please 
refer Ketchum [3] ). 

The modified Thomas algorithm v^Bs employed to solve 
the system, and the various steps of which are presented in 
the Chapter 2, 

By exploiting the sparsity a reduction in the 
computations by a ratio of 3,013 was obtained, as only 
62961 operations were performed instead of^ 189702 standard 
operations. The comparison and details of the standard 
and improved operation counts have been presented in the 
Table 4, 

There are 46 upper diagonal, 53 off-tridiagonal, and 
hence a total of 99 P-blocks in the matrix, after the forward 
substitution. Each P-block is of the order 9x9, As 
discussed earlier only a 9x5 matrix for every P-block 
needs storing, since the remaining elements are zero, and 
hence a 44% saving in the storage may be realizedi 



10 PLATES 
P= 20 atm 



10 PLATES 
P= 20 atm 


15 PLATES 
P= 4 atm 


12 PLATES 
P= 15 atm 


FIG 5 THE ARRANGEMENT OF THE INTERLINKED-CDLUMNS 
(PROBLEM 4) 
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A? Eb C? 
A3 B3 
A4 


C 


C3 

a C4 

Afi Bfe 

A? 


Ce 

C 7 

Ae 8 b Ce 
Ag Bfe C 9 

Ajo Bic 


Bjj C]j 

Ai 2 Bj? C ]2 
Ai3 Bi3 Cj3 
Ai 4 Bi 4 Ci 4 

Ai5 Bjs C35 
Ajfi Bj 6 C16 
Ai7 B]7 Ci7 


C n.ii 


Asa,! 


Aib Bj 8 CiB 
Ai 3 Bjs Cig 

Ago ^ 


621 C 21 
A22 822 C22 
Az3 623 

A24 024 C24 
A25Q25C25 

A26Efe6C26 
Ae? ^ C27 
A28 ^ Cze 

A29B2S C29 
Ago Bai Cao 
A 31 ^1 C 31 
A32B32C32 

AsaBsaCas 
A 34 B 34 C 34 
A 35 B 35 


A 402 ] 

A 41 ^ 


Bae C36 
A37 B37 C37 

B38 Cse 
A39 B39 C39 
A40 Bto C40 
A41 B41 
A42 


C41 

B42 Zi2 
A43 B43 C43 
A44 B44 
A45 


C *4 

B45 C45 
A 46 B46 C46 

A47 B47 


FIG 6 the structure OF THE JACOBIAN MATRIX 
FOR PROBLEM 4 
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TABLE 4 

Comparison of Operation Counts for Problem 4 
Number of components = 4 
Total number of stages- 47 


Operations 

Number of 


Standard 

Improved 


times 


operati on 

operation 


performed 


count 

count 

Op 1 

42 

x729 

= 30618 

x225 = 9450 

Op 2 

2 

x729 

= 1458 

X 81 = 162 

Op 3 

43 

x729 

= 31347 

X 41 = 1763 

Op 4 

38 

x729 

= 277 02 

X 41 =1 558 

Op 6 

2 

x729 

= 1458 

X125 = 250 

Op 8 

36 

X729 

= 26244 

x225 = 8100 

Op 1 0 

1 

x729 

= 729 

x405 = 405 

Op 1 1 

25 

x729 

= 18225 

x225 = 5625 

Op 12 

1 

x729 

= 729 

x225 = 225 

Op 1 5 

45 

X 81 

= 3645 

X 9 = 405 

Op 16 

2 

X 81 

= 162 

X 25 = 50 

Op 17 

66 

X 81 

= 5346 

X 45 = 2970 

Op 1 8 

40 

X 81 

= 3240 

X 45 = 1 800 

Op 19 

47 

X 81 

= 3807 

X 81 = 3807 

Op 20 

47 

x729 

= 34263 

x546 = 25662 

Standard 

matrix 

multiplica- 

tion 

1 

X729 

= 729 

x729 = 729 

Total 



1 89702 

62961 


S.Q^C. 

I.o.c. 

= 3 

.01 3 

, 7 . ■ 
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It may be noted here that out of the 99 P-blocks, 

15 are null matrices Pis^20* ^21,47* ^22,47* 

and ^25^47 through ^j) and hence need not be stored, 
thereby leaving only 85 P~blocks v\tiich must be stored. 
Therefore, only 3825 (= 85x9x5) elements are stored instead 
of 8019 (= 99x9x9), and a 52% of actual saving in the storage 
is obtained. 

Problem 5 ; Consider the system of interlinked columns for 
separating a mixture of 3~components using 3 columns which 
are ccnnected as shown in the Figure These columns are 
assumed to be split into section 2a having 4 plates, 2b 
having 3 plates, 3a having 4 plates, 1 having 4 plates, 
and 3b having 4 plates, and are rearranged in that order. 

The Jaccbian matrix for this system is presented in the 
Figure 8, The modified Thomas algorithm was applied to 
obtain the solution. 

In this case, the sparsity exploitation has resulted 
in a reduction in the computations by a ratio of 2,744, 
as only 13393 operations were performed instead of the 
36750, The comparison and details of the operation counts 
are presented in the Table 5, 

There are 18 upper diagonal, 11 off-tri diagonal, and 
hence a total of 29 P-blocks, each of order 7x7, in the 
Jacc±)ian matrix, (P| through P^ ^4 1 2 P4^-|o» 

Pl1 16 P -|4 ^ 5 ^ shown in the Figure 5, Since only 




FIG 7 AN ARRANGEMENT OF INTERLINKED COLUMNS 
(PROBLEM 5 AND 6) 
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8 THE STRUCTURE OF THE JACOBIAN MATRIX 
(PROBLEM 5 AND 6) 



TABLE 5 

Comparisai of Operation Counts for Problem 5 


Number of components = 3 

Total Number of stages = 19 


Operation 

Number of 

times 

performed 

Standard 

operation 

count 

Improved 

operation 

count 

Op 1 

18 

x343 = 6174 

X112 

= 201 6 

Op 3 

30 

X343 = 1 0290 

x25 

= 750 

Op 7 

10 

X343 = 3430 

x112 

= 1120 

Op 9 

1 

x343 = 343 

x196 

= 196 

Op 14 

18 

x343 = 6174 

x112 

= 2016 

Op 15 

19 

X 49 = 931 

X 7 

= 133 

Op 17 

40 

X 49 = 1960 

x 28 

= 1120 

Op 19 

19 

X 49 = 931 

X 49 

= 931 

Op 20 

19 

X343 = 6517 

x269 

= 5111 

Total 


36750 


13393 


75 


a total of 812 (= 29 X 7 x 4) elements are stored instead 
of the 1421 (= 29x7x7), a saving of about 43?o in the 
storage requirement is obtained. 

Problem 6 ; The same system of the interlinked 
columns described in the problem 2 vas solved mth a 
6~component mixture. The structure of the Jacobian and 
the method of solution remains the same, but the order of 
each submatrix in the Jacobian is 13x13 now; 

In this case, the sparsity exploitation has resulted 
in a reduction in the computations by a ratio of 3*21 , which 
is more significant than that in the problem 2; Only 69811 
operations were performed instead of the 224094, and the 
details are presented in the Table 6* 

The total number of P-blocks is 29 (same as in 
problem 2), A total of only 2693 (= 29x13x7) elements are 
stored instead of 4901 (= 29x13x7), and thereby a saving of 
A-6% can be realized in the storage requirements.* 
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TABLE 6 

Comparison of the Operation Counts for Problem 6 

Number of componoits = 6 
Total number of stages =19 


Operation Number of Standard Improved 

times operation operation 

performed count count 


Op 1 

18 

X2197 = 39546 

x637 = 11466 

Op 3 

30 

X2197 = 65910 

X 85 = 2550 

Op 7 

10 

X2197 = 21970 

x637 = 6370 

Op 9 

1 

x21 97 = 21 97 

x1 1 83= 11 83 

Op 14 

18 

X2197 = 39546 

x637 = 11466 

Op 1 5 

19 

X 169 = 3211 

x 13 = 247 

Op 17 

40 

X 169 = 6760 

X 91 = 3640 

Op 19 

19 

X 169 = 3211 

X169 = 3211 

Op 20 

19 

X2197 = 41743 

x1562= 29678 



224094 

69811 

I.O.C. ‘ 
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CHAPTER 7 
CONCLUSIONS 


In the algorithm presented in this work the sparsity 
and the structure of the submatrices was exploited in the 
various matrix multiplications and inversions, and thereby 
a significant reduction in the operation-count and storage 
requirements was obtained* The reduction in the computations 
becomes increasingly significant with the increase in the 
number of components. 

Theoretically, a maximum of 3,43 (= 24/7) times 
reduction in the computations may be realized in the single 
column problems if the number of components is very large* 

In general, a reduction in storage requirement by a ratio 
of (c+1)/(2c+l) was also realized. The exploitation of 
sparsity in solving the single column problems by the 
conventional block Thomas algorithm has shown 2,676 and 
3,02 times improvement in the operation-count, in the test 
problems 1 and 2, 

An efficient approach to solve the problems with 
the intermediate tray specifications has been proposed and 
shown to be more advantageous than the one proposed by 
Hofeling and Seader [2] , in both the computational and 
the storage aspects. For instance, in the problem 3, the 
proposed approach has been proved to be 3,877 times 
computationally efficient and 64^^ more storage-saving than 
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the Hof eling-Seader approach. 

In the solution of a system of interlinked columns 
by the modified Thomas algorithm, the exploitation of 
sparsity has resulted in a 3,013, 2,744 and 3i21 times 
improvement in the operation count, in the problems 4,5 
and 6 respectively. Though this ratio varies from problem to 
problem, and improves with the increase in the number of 
components, in general it is close to 3, 

The saving, at least by a ratio of (c+l)/(2c+l) 
in the storage requirement is also realized in the systems 
of interlinked columns, vhlch significantly improves 
further if some of the of f-tri diagonal P-blocks are null 
matrices. 

The proposed method of the -sparsity exploitation 
can bo applied to the other related techniques, for instance, 
the convergence domain extension methods proposed by 
Vickery and Taylor [8] ; and the other methods of solution 
proposed by Kubicek {7] , and Stadtherr [ 5 ] •; 
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VAPOR ENTHALPIES POP TTIE TEnPERATURS RANGE OF - 25®F To 40°P AT P « 800 Ib/in abs. 
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APPENDIX B2 


SA-MHiE H^OBLEM No. 2 
DAm 

K values 


Ma t er ia 1 Temperature 

100°P 200°F 


A 

500.0 

550.0. 

B 

1 .50 

1.8 

C 

0.90 

1,00 

D 

1 .0 X 

10“^ 1.5 X 10"^ 


Molar liquid 

_ 

enthalpies, 10 cal/mole. 

A 

0.01 

0.013 

B 

0.50 

0.33 

C 

0.40 

0.44 

D 

1 .50 

1.90 


Molar vapor 

enthalpies, 10 cal/mole 

A 

1.00 

1.002 

B 

1 .80 

1 ,82 

G 

2.00 

2.03 

D 

5.75 

5.95 



; lifKKiilNKKlU COhUms 








PROGRAM INTGNK 
INTEGER CfCPl,T«OC,CT 

DIMENSION.,JCO«D(23;iT5fPEC(2)taReBGC2)|inPER(2) 


I Ml «» 1 


l IXt6),llfC6) ,GORVce»)rRATIOCi),aNOSVcS), JNOSliCaS, 

I aVAPrCl)/TFV(l),pEE6vU»I.)f4L.lQFClI .TFGCl) fFEE^IjCl,!) , 
i JSPC(4),£Bt(4).SpC|^AL(45,l6aMF(4h ■ ' 

1 QQap.SVtl2,25»St(12f2)eSSVCl2),SSG(l2)#T(U),ETACl2) 

5 gJUIS • 

DIMERSIufi ACg,9),Bcy,95,CCC9,9),A0fF(9,9),C0FF(9,9) 
DIMENSION A£.FAcy;&;,6ETA2C9,5I,eETAlC95,VEClC9J 
DIMENSION FC47;95#H8!>, 9,5), ac47^9) ' 

t WIM « IM IM IM in «•» fM -Ml im «i» IM w iM* ^ Mt A * M IM' Ml fM IM IM «i» i» P» 

0PENCUNlT=3J,DEVICEs'DSM',FIDE«*'INP,iN') 

OPEN t:UNITs51,DEVl€Es*'DSM*,FItEs'TH, OUT') 

REAOC33,*)N 

REAOC33,*JC,GPl,T«OC,CT 
REA0C33,*)NOFF-NAliOVE 

REAOt33#>^)NCOND,NREBD,NUy,NOL,NVAPF,NLIQF,NSPC 
IF(NC0 n 6,GT,1)G0T0 810 
DO -800 Isl-MCOND 

READ(33,?)5coNDCI),ITIPCC{I) 
lF(«REBt.,Dr,i)GaTO 811 
DO -801 lal-WREBD 
REAOC33,#)5HEBDCI),irypERCl) 

IF(NSPCrDT, 1)60X0 812 
00-802 Is1,NSPC 

8EA0C33,»)JSPC(I),DBD(I),SPCVAD(I) ,1C0MP£1) . 
IFlNi1V,tiT.l)G0X0«l3 ' 

REA0t3§,*)taN0SyCl) ,I=lrN0V) 

IFINOD.DT.DGOTO -814 
READ133,»)CDN0SDC1) fI*l,N0B) 

IF t NOFF . LX ; 1 ) 6UX0 815 

00 803 isl.NOFF . , , . 

REAO£33,»JIX(X)ai(I),LURV(I),RATlOCl) 
IF£NVAPF,LT,1)G0TU 816 
DO 804 I=1,NVAPF 

BE AD C 3 3 , » ) J V APF 1 1 ) , XF V 1 1 ) , (FEED V ( I , KK ) , KK*! , C) 
IFtNLlQF.LX.DGOIO 817- - 

DO '805 I=1,NLIQF 

READ C 3 3 , >^ ) aLIQF C I ) , TFL 1 1 ) , C FEEDH I , KK ) , KK^l , C ) 

CONTINUE 

REAOC33,*)(QQCJ),asl,N) 

READC33,*)CSSVCa),J=l,N) 

READC33,»)CSSLC0),J*1,N> 

REAOC33,>l')CETA£J),J = l,N) 

KEADC33,=l')CfCJ),a»l,N) ' ^ 

READC33,X‘)C(SV£a,l) ,iai,C) ,0=1, N) 
HEADC33,*)£(5LCa,l),I=l,C),0=l,N3 

I iP in IM »M IM Ml Ml ■»# IM »M *» i» Ml .in w ^ I* 

XCHK=0 
KCONDSO 
KREBL*© 

NTOfsNCOND+NREBL 
DO IID KK=1,NSPC ^ ^ , 

IF£LBL(RK),LT,3)G0T0 IIU 
JsJSPCCMK) 

DO 311 MM«1,NC0ND 
IFCvI,EQ.JCDND(«M))GOTO 113 

111 CONTINUE 
DO 112 «M=i,NREHL 
IFCJ«E0.JKEBL(MM))G0T0 114 

112 CONTINUE 
GOTO 115 

II J KC0ND=KC0ND+1 


800 

810 

801 

Ul 

802 

812 

813 

814 

803 

815 

804 

816 

805 
81 1 



2 


114 

lib 

110 


lib 


117 

118 
yui 
902 
908 
904 
119 

I ifmmm 

145 


14? 

146 


151 


152 


158 


154 


GHTfl 115 

KREBt,»KKEHL+l 

1CHK=J,CHK+1 

CONTINUE 

1SXP=0 

XFCKCONU.EQ,NCOND3GOrO 116 
IS’XP=1 

XYPE 901,KCOflD,NCONO 
iF(KRg;Bo.fc;Q,f<Ri;»i.)Goro u? 
ISTP=1 

XyPE 902,KREBf4,r4RKBl. 
IF(,ICHK,EU.«fOf IGUTO 119 
XFCICHK,L'l',«XOT3GUTU 118 
!l*YPE-9a4,lCHK,NT0X 
STOP 


XYPE 90i,lCHK,Nt0T 
STOP 


FORMAT (I 
1 HAS 

format t I 

1 HAS 

FO^MATtl 
1 BUT 

FORMAT (1 
1 (' BUT 
IFtlSTP. 


X,'-— EKRUR— VlX.'IffFORMATION ABOUT *’,12,' CONDENSERS' 
' BEEN PKOVIUEO whereas TMERE'ARE', 12,' CONDENSERS') 
X,'— error— •'/IX, 'INFORMATION ABOUT ',12,' REBOIbERS' 
BEEN PROVXOEU WHEREAS MERE ARE ',12,' REBOiLERS'J 
X,'-— ERROR— '/IX- 'NO, OF SPECIFICATIONS SHOULD BE ',12 
',12,' rave ttEEN PROVIOED'/tXi'PROBtiBR UNDERSPECIFIED') 


',12,' rave ttEEN 
X,*»— error— '/IX 
',12,' HAVE BEEN 
EOU)STOP 


ROK— '/IX, 'NO, OF SPECIFICATIONS SHOULD BE '.I 

HAVE BEEN PROViDED'/tXi 'PROBLEM OVERSPECIFIED*) 

£!> --■■■ ... - .« ..... .. .... 


DO 145 iJ(sl,N 
RD(J)=1, 

RVCJ)*!. 

DO *146 KKsi,KOrr 
JalYtKK) 

IFCDORVCKKJ,EO,0)GOTO 14? 

RVta)sRV{J)-RATIOtKK) 

GDT0146- 

RL (;J)-RD(J) -RATIO! KK) 
CONTINUE - 
00 151 J=1,N 
8VtJ)=0. 

BUU)=D* 

DO 151 1=1, C 
8V(Jj=BV(j5+SVCa,l) 

6UCJ)=fiDCU)+5ECU,l5 

DD-200 U=1,N 

jMi=a-i 

jpi=j+i 

TvJ=TtU) 

SDO="Ci ,+850(1)) ) 
SD1=-(1,+SSV(U) ) 

DO 152 KK=1,NC0N0 
IF(J.EQ.aCONDCKK))GOra 180 
CONTINUE 

DO 158 KKsi,«REBL 
IF(J;EQ,UREBDCKR))GaTO 190 
CONTINUE 
CADE ENVCTJrHV) 

CALD ENDCTU,HU) 

CAJbli FINDKCTU,AK) 

SUMVaO. 

SOMLaO. 
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DO 15)5 KK=l,NOL 
IF(J,EQ,JNOSl,CKK) JGOrO 170 
CONTINUE 

CADIi ENUCTta-l) ,HD) 

SlJHML=0.0 

DO '171 K,KsjfKaV 
IFlJieQ.UNOSV(KK)lGOrO ;«oo 
CONTINUE 

CALL EN(r(T(J+l) ,HV) 

SUMHV= 0 , 

DO 173 Isl.C 

rcj,i)aF(u:i)+sv(api.i)»RVcapi) 
SlIMHV=SUMHHHVCI)»SV£aPltIJ 
FCU,CT)=Fia,CT)-»-SUHH\f*RVtapi) 
GOTO 200 


I ifMlW W|li||ill«kll» 1 MWll» 

t CONDENSER PART 


DO 181 Isl,C 

rcj.i)=svijPi,i)*R 7 caPii+s 0 i»svcj,i)+SDO*SLca,i) 
irClTYPECtKR),EQ,0)GOTO 182 
DO -184 1=1, C 

FCa»l+C 3 sBVtJ)»SliiafI)/BLCJ)»SVlJ,l) 

GOTO 200 ' ■ 

CALL FINDI((TJ,AK) 

C 0 EFlsETA(D)»B¥(ai/BLt 31 

C0EF2=(l,-EfACJ))»BVCU)/&VcaPl) 

DO 183 1=1, C - . 

FC U , 1+C J sCOEFl ♦AK 1 1 ) l^SL t J , I ) -SV C J, 1 5 +CQEr2»S¥ C DPI , I) 
GOTO 200 

p»«p mmm wn «ii mm mi w . • 


t REBOILER PART 


DO 191 1=1, C 

FCJ,l)=SLta«l,I)*RL(aMl)+SDl»SV(J,I3+SD04SL(J,I) 
ir(ITYPERtKK3,EQ,0)GaTO 192 ' ' 

DO 194 1=1, C 

FCJ,l+C)=BV(.J)»SLtJ,lJ/»L(a)-SV{J,I3 
GOTO 200 

CALL FINDACTJVAK) 

COEF=8Vta3/BL(a) 

DO 193 1 = 1 , C 

FCJ,l+C)=c 6 EF»AKCl)*SHJ,I 3 -SVCa,l) 

CONTINUE * - 


S FEED PART OF THE DISCREPENCY FUNCTIONS 


DO 201 KK«1,NL1QF 
asJLIQFlKK) 

TF=TFLCRK) 

CALL FiLCTF»HL) 

DO 201 1=1, C 

F(a,I3=rtD,I)+rEEDLtKK,13 
FCJ,CT)=ria,CT)+HLCI)*FEE0L<KK,l3 
DO 202 KK=1,NVAPF 
a=JVAPFlKK) 

TFsTFVCAKJ 
CALL ENVCTF,HV) 

00 202 1=1, C 

FCJ,l)=FtDil} + rEED¥tKK,l) , 
FCa,CT)srta,CT)+HVa)»FEEDVCKK,I) 

» » III* »«* «» ,1# iw w «» <*' ')|i» fNi |i» «• 'fig «i* i# 


: CONTRIBUTION OF INTERLINKS TO fHE OISCREPENCY FUNCTIONS 
DO 210 KR=l,NOFr 



EQ,0)(J0fa 

HVJ 


211 


21J 

211 

212 

2ia 

! M iw iit iMi int M 

1 EASf 

Jr mmmmimm 


JK=II(KK) 

I'KsTCJK) 

IFCIitlRVlKK), 

CAlilj ENVCJ'Kr 
SUMVsO, ■ 

DO 213 1=1, C 

SUMV=SUMV+liV 
F(d,eT)=FtJ, 

GOTO 21M 
CALL ENLCTK, 

SUMLSO, 

DO 212 I«1 ,.C 

f ^ 1 ^ A P t ^ ^ ^ 0 c K K 3 

SUML=SUML+SLt«JK,I3»MLtI3 

f;Ca^CT)=Fld,CT)+SUML»RATIOCKK) 


3+SVtaK,l)«KAT10(KK) 

(I)>»'SV(aK,13 

CT)+SUMV»WATiOCKK) 

HL3 


K H* f* •• «•«■■■» HI ■ 


220 

221 

222 

223 

225 

227 


225 

905 


230 


231 


233 

232 


235 

237 


TERM OF THE EQOILieRlUM OISCREPEMCY FUNCTION 

DO 240 asl.N 
XF(ETACd),£Q,l,3G0T0 240 
DO -220 KK»1,NC5 nD 
ir( J,EQ,aCDNDCRA) JGUTU 240 
CONTINUE 

00 221 KK=I.NREBL 
IF(d,EO.JKEOLCKK:)3GOrO 240 
CONTINUE 

KK=l 

IF(d.NE,IXCKK>)GOfO 223 
IFCLORVCKKI.Ea.lJGOTO 230 
KKsKK+l 

IFCKK.LE-NOFDGOTO 222 
00-225 RK*1» NOV 
IF(U,EQ.dNDSVCKK)3GOTO 226 
CONTINUE 

COEFFs ( 1,-ETA CJ3 3 »BVtJ)/'BVCd + 13 
DO 227 lsl,C 

F(J,l+C3=FCd,I+C3+COErF»SVCa*H,I) 

DEROMt J3 = KVt*J+i3=rBV(,J+t3 - 
GOTO 240 
TYPE 905, J 

FOKMATUX,"— error— '/ iX, 'STAGE ',12,' ODES NOT HAVE ANY *, 

1 'VAPOR INPUT ANO HAS WOT BEEN OEFINE0 AS A REBOILER') 

STOP ■ - 

LLsKK+1 

JKsIYCKK) 

DENOMCJ)=BVCJK)*RMTlDtKK) 

CDEFFsCl,-ETACd))»BVCJ) 

DO 231 1=1 ,C - 

VECT0RtI)=SVCJKrl3»RATI0(KK) 

IFCLL.GT.NOFFJGOrO 232 
DO -232 KK=LL,NOrF - 
IFCJ,NE,iXCKK))GOTO 232 
iFCL0RVtKK),E0,0)G0T0 232 
K "** C K K 1 

DENdMCd)=OEMOM(J)+BVCdKJ»RATIOCRK) 

DO 233 1=1 ,C ' ■ ■ ■ ■ 

VECTORCl)=VECTOR(l)+SVCaK,I)»RATIOCKK) 

CONTINUE - • - - 

DO 235 KK«1,N0V 
IFCU.EQ,JNOSVCKK)3GOTO 238 
CONTINUE , . ^ 

OENOMCd)=DENOMCJ)+BV(J + l)»RVCd+U 
00' 237 C " * ' . - . . 

F(J,l+C)=F|o,I+C)+COErF*C^ECrORtI3+SVCa+l,I3»RVCaTl>)/OENOflCa) 
GOTO "240 ' ■ ' , 
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IB DO 2iy 1 = 1, C 

'r !! r.«r!!r"* 

SPECIFICATIONS PAKf 

» !P» *» «» «» IPB ((I* ^ » «P *» IW III, ,«P ^ ^ ^ ^ IP! ^ ^ ^ ^ ^ ^ ^ ^ ^ 

K K 5» 3 , 

'0 J=dSPC(NK) 

api=a+i 

TJaTCJ) 

SD0=«(l,+SSl<(J)3 
SD1=«(1,+SSV(j5 ) 
lFC0BL(KK),l4f ,&3G0I0 260 
DO-251 MM=1,NC0ND 
IFtJ;EQ,JCOND(MM))GOTO 255 
>1 CONTINUE 

00 252 MM= 1 ,NRE 6 D 
iF(a;EQ,dRE 6 LCMM))GOTO 253 
>2 CONTINUE 

GOTO 260 
>3 QQ(J)=0, 

CAbU BNIU(T(a-13 ,Hl4) 

CAliD ENVCTJ.HV) 

DO 254 1=1, C 

>4 QQ Cd)=QQ(u5 + (Nil (!?•«¥( I) )4SV(d, 13 

QQ(J)aQa(U)4SDl 
GOTO 260' 

>5 QO(J)=0,* 

CAbL ENV(T(U+ll,HV) 

CAIiD ENli(Td,Hli) 

IFCITyPEeiMM),E«,0)GOTO 256 
DO '258 1=1, C 

>8 Q 0 Cd)= 0 Q(a 5 +(HG(Il-HV(IJ)»SVCaPl,I) 

QQCJ3=Gl«Cd5»RVed-»-i) 

GOTO 260 - 

>6 DO 257 1=1 ,C 

37 QQld)=QO(d5+CHV(IJ*HbCIJ)»SL(d,13 

QQCU)s 00 W)»SDO ■ 

&0 GOTO (10, 20, 30 d 40 #50, 65, 60, 65, 70, 75, 80, 85, 90) ueiiCKK) 


260 

^ mmm* 

t NO 
10 




IS NOT A CONDENSER’ 


CONDENSER SFECIFICATIONS 
DO 11 MMsI.MCOND 
irCU,EQ,JCONO(MM))GOTO 12 
CONTINUE 
TYPE 906-0 

FORMATCIX,'— error— ■•'’/ iX, 'STAGE ',12,' IS NOT A CONE 

STOP 

QQ(d)=0. 

FCd,CT)a0, 

CAliD ENli(TJ,Hli) 

CAG6 ENVCT(0+1),HV) 
irCIfTPEClNM),EQ,U5u0T0 14 
D0'13 I=i,C 

QOCd)=QaCU3+SV(OPl,l)4(HLCI3*HVlI)) 

r(3,cT)=r(j,cT3+Hb(i)»cH¥CaPi)*svcjPi,i)+SD04si,ca,i3 

1 +5D1»5VIJ,1)) , " 

aQ(j)=QO(d)’f‘RVcapi) 

GOTO 265 • 

DO 15 1=1, c. 

Jtc:*r)+Hvci3*cKicdPt)^l{caPi,i)+SD0»5ii(iJ,i)) 

QQ(J3=QaCJ)’SD0 
CAbU ENV<TJ,HV) 

SUM=0. ' 

DO 16 1=1, C 
SUMaSUM4HV(i34SVCa.I) 

FCJ,CT)=r(d,CT)+501*'‘SUM 



6 




GOTO 265 


t NO RtbOILER specifications 

20 DO 21 MMsl.WREBO 
irCJ,Ea,aKEBE(MM)lGOT0 22 

21 COJITINoe 

rypE 90 /.a 

907 FORMATCXXr'— ERROR— Vix, 'STAGE ',12, 

STOP ' - ■ ■ 

22 00(03*0, 

FCJ,CT}={). 

CAOli ENli(T(J-13 ,H0} 


IS NOT A REBOIIiER'3 


26 


CALC ENVCT0,MV3 
IF(ITypER(HH3*eO,03GOTO 23 

00"2lti 

OQ(03=0a(0)+CHECI3-H¥(I3)*SVCa,l3 

f ^ ^ « ST W ' ^ ^ ^ ^ ^ 1 ^ ^ C ^ > ’♦'SO 1 


P vw Jl-f '.V 

1 ♦SVCJ,I3) 

00(03=501*00(03 
GOTO 265 

23 DO 24 1=1, C 

aQC03*QO(03+(HD(I3*HVtI33»5V(a,I3 

24 F(O,CT3*F(O,CT3+HD(13»CKD(OH13*SD(OMl,13+S01*SVCa,I)3 

0Q(03=SDl*QO(03 - 

CAUL END(T0,HD3 
5UH=0, 

25 i8MiioH*kfl3»SL(O.I3 

r CO^C||=F (0 , CT3 +SDO*SUM 

• «V« IK W lip MN«p«»ipW p« IP MR 

{ CONUENSER-HEAT-DUTT SPECIFIED 

30 DO 31 MMalrNCOND * 

IF(O.EQ,OC6NDCMA3)GOrO 32 

31 CONTINUE 
TYPE 90«,0 

908 FORMATdX,'— -ERROR— '/IX, 'STAGE ',12,' IS NOT A CONDENSER') 

32 FCO,CT3=SPCVALtKK3 
QQ(03=SPCVA)U(KK) 

CADL ENVCTCO + DfHV) 

CALL ENL(TO,HL) 

SU«HV=0, 

DO 33 1=1, C 

33 SUMHV=SUMHV+HVCI)*SVCaPl,l) 

F(J,CT3=F(J,CT3+RV(aP13*SUMHV 


34 

35 


r vw «v.A j-f vufv.^, 

CALL ENV(T0,HV3 
IF{ITYPECtMH3,EQ,03GOTO 35 
DO ■ 34 1=1 

rCaiCT3=F(0,CT3+HL(l3*(5D0*SL(0,I3+SDl*S¥(0,I33 
GOTO 265 


SOML=0, 
simv=o, 

DO 36 1=1, C 

stmv=suMv+Hvci3*svca,i) 

36 SUML=SUML+HL(I3*SLC0,I) 

r(0,CT)=Ft0,CT3+SD0*SUMLFSDl*SUMV 
^ GOTO , , , , , , 

j"REioiLiR"HEAT-DUTY'"sPEGiFi:iD"' 

40 DO 41 MMsl-NREBL , 

IF(O.EQ,OR£BLCMM33GOTO 42 

41 CONTINUE 

TYPE' J 

909 FORMATUic,'— ERROR— '/IX, 'STAGE ',12,' IS NOT A REBOILER'3 

■ ' STOP ’ ' ' ■ *' 

42 F(0,CTpSPCVAL(KK3 
QO(03=SPC¥AL(KK3 
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4i 

44 

45 

46 


CALL. ENL(T(J-1) ,HL) 

CALL ENVCTJ^HV) 

sim=o. 
on 4i 1=1, c 

SUM=&UM4HV(l)fSV(a,l) 

FCU,CT)sFCJ,C’r)+50i*SUM 

yQ:44_i=i;c _ 


PC^gCT)=F?0,CT)+HLa)'^(«LCaMUfSL(0«l,IHSD0*SL£a,I3 3 


StJMsO, 


SUMlaO, 

CALL ENL(TJ,HV) 
DO 4b 1=1, C 


SUM=SUM>«'HVCl)»SL(a,I) 

&l.lMl=SUMl+HLCl)»SLCaf«l,l) 

FCJgCT3=FlJ,CT)+5O0*SUM<-RLcaMl)4S0Ml 


t REFLUX RATIO SPECAFIEO 

50 DO 51 MMsl.NCORD 
IFCJ.EQ.JCONOeMM} JGOrO 52 

51 CORTINUE 
TYPE 910,0 

910 FORMATdX,'— error— V lX, 'STAGE %I2,' IS NOT A CONDENSER*') 

52 FCa,CT)»BLCO)-BV(a)*SPCVALCKK) 

GOTO 265 ' 

t KEBOILER RATIO LV/B) SPECIFIED 

55 DO 56 NMssI.nrEBL 
IF(D,EQ,aREBLCMN))GOTO 57 

56 CONTINUE 
TYPE 9ll,«J 

911 FDRMATCIX,'— ERROR"— VIX, 'STAGE ',12,' IS NOT A REBOILER') 

STOP" ' "■ 

57 FCa,CT3*BV(03"BLCU)»SPCVAL(KK) 

GOTO 265 

i TEMPERATURE Ur THE STAGE AS SPECIFIED 
BD F(d,CT3=TCa)-SPCVALtKK3 

GOTO 265 ■ 

^ Mi MN KM Ml Mi Ml Ml Mi KM Mi Ml VM Mt Mi 'Mi MW Ml m MM m, Mil Mi M<t MW MW Ml MW Ml MW Mr MW iM WM mu Ml tM MW WM MW iM MR WM WM Mi WM Ml 'MW RM Ml WM WM MW MW .m|. IM WM Ml WW 

S TOTAL VAPUR-FLON VCJ) LEAVING STAGE J IS SPECIFIED 
t t OK "D« IK CASE OF A CONDENSER 3 

Ml «M Mi MR <M IM'«MM»4MMWMRiMMBiMMWMiiMMlWMiM'M»M»fMiMMifmiMWM'iMiMiMMilWMlMttMRMWMM''MiMWfM M»«W'MRM»W^ 

65 FtJ,CT3SBV(0)-SPCVAHKR) 

GOTO -265 ■ ■■ ■ ■ . - - 

f «MWM'tlWllM«M<MiM«M'MWMMMW'MWMW|MWMMWWMMWMM|MMRMW*MMWMip|lMiiMifMM'M*MK<MMW'WMMRWMpMMWiMMiMI*MMiMM 

S TOTAL LIQUID-FLUN LtJ) LEAVING STAGE J IS SPECIFIED 
t L OR "B* IN CASE OF A REBOILEK 3 

f Mi i# Ml IM MW IM IM » MW •• Mw Mi MW IM M» WM «M «M ^ MW M» MR IM WM MW M# «» Ml MW WW-' MW Ml '*M Ml ^ ,MW ,M» WM MM Ml W MW Mi .WM «M MW Ml IM Ml M» *» • MW Ml IM MW MW MW . 


70 F(a,CT)=BL(0)-6PCVAL(KKJ 

GOTO" 265 

:"cOMPONiNT-FLOiiIyAPOR 3 *vl 1,13*0? CoSpONENT I LEAVING STAGE 3 SPECIFIED 
t C OR DCl) IN CASE OF CONDENSER 3 


75 


1I=IC0HP(RK) 

F(0,CT)=SVCU,II3-SPCVAHKK3 
GOTO 265 - 

ilM'M»M»«MiMl|WiMWMMW:WMMiMWM»MWWlW'MIMtWM'iMMW'|MMI'MiMWwM'M»'MW 


S CUWPUNEKT-FL0WtLI«UID3 LC3,13 OF COWPONENT 1 LEAVING STAGE 3 SPECIFIED 
S E OR BC3L3 IN CASE OF RE601LER ) 

80 IlsICOMPtEK3 

FC0,CT3=SLCD,II3-SPCVALEKK3 
GOTO' 265 ' 

i VAPDR-MOLE-FHACTION ytS,!) OF COMPONENT 1 LBAVIMG STAGE 3 SPECIFIED 




OR XdCi) IN CASK UF CUNOENSER 3 
1I^^C0«PC^^J 

^fJgCT^=SVCJ,IIJ-SPCVASlKK3*BV£03 

• «» «t «»*» w# »*,«« Iiw ^ ^ <|I^ w ^ flit ■•» iw. It B* ^ «»*»» ,f^i » ^ ^lai w <»« »pt- • w ip »»•■<» lii iw «|^ Ii» 'ap 1^ 

1 leaving stage i specified 

IIsICOMRCKKI 

FCO,CT3=SSCO,II)-SPCyAl,UK3^'aLCvJ) 

P ^ ISI» ^ p» p» |« IIP* l|(p> i» i*^ 111 ^ ilil ,1(1 p» ^ Ml* p, ip. ^ ^ ^ |g> P, ^ I*. ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

1 KKsKK+J 

1F(KK.LE,N5RC3G0TU 270 

cntiTiNue 


DO 400 Usl.N 

«(RIfeC51p*5o,CFCJ»IJ#I-t,CT) 
wRiiEC5i,»3c6Qca)fa»ifNj' 
iKRITE(52,»3CRVCa)pJ»l,N) 
WRITEC52,»)tRl»CJ)^ J=i»NJ 
WRlTEC52,»)CDEWUMCa}pJ=l,N3 

««|l»«Mil»f|tt|i»ii|»M4W«MlRVMll«»«»WM»«»-|W«i»lWMMkW 

IZl=47 

IZ2S57 

BliD»2O0« 

RL(13=200,/Bl*(20) 

RLC413sl,''RD(13 


SOD 


901 

1 mm 


902 


Jsl 

CALO 

CALL 

CALL 

CALL 


TRIBCJ) 

TRICtJl 

OFrctJ,2O,T(2033 

INVPRICB) 


CALL BCMULC4) 

CALL erMULCJ) 

CALL 8A«ULC47,B,COFr3 


DO 501 Oa=2»9 
OaJd 

jz=izi+oj-i 

CALL TRIAtJ) 

CALL TRIBtJJ 
CALL TRiCtJl 
CALL APVMOLta-1) 

CALL AQMULCO-1) 

DO 500 1=1 pCf 
rCJ£l)s:FCD,l)-VEClCl) 

DO 500 Ksl,CPl 
eH,K+C3=BU,K+CJ*'ALFA(l,K3 
CALL INVPKTtS) 

GALL BCMULCO) 

CALL BFMULCO) 

CALL APLMULljZ-1) 

CALL MBAPLCaZl 
CONTINUE ' 


JSIO 

JZ=S6 

CALL ABSBOT 
CALL TRIAIJ3 
CALL TRIBtJ) 

CALL APVMOLU-iJ 
CALL AOMULCa-n 

VEClCl) 

■l»C3-ALFACArK) 


DO 502 1=1 ;CT 
F(J,I)=FCJ, 13- 
DO 502 Ksljcpi 
BCI.K+CJbBU.K 
CALL BFMULC^J 



9 


50fo 

504 

J rnmm 

507 

511 

509 

510 

508 

512 


CALL APLMOlCJZ-1) 
CALL MBAPL(02) 


J=ll 

JZ=57 

CALL ABSTOp 
CALL TRIBIJ) 
CALL TRICtJ) 


CALL 

CALL 

CALL 

CALL 

CALL 


INVPKf CB) 

BCMULf J) 

BrMOLfJ) 

gFrAtJ,47,T(47)) 

8AMULC57,B,C0Fr) 


DO 504 00*12,19 

JsJJ 

JZaJZ+l 

CALL TRiACJ) 

CALL TRIBCJJ 
CALL TRICCJ) 

CALL APVMULta-1) 

CALL AQMUL(iJ-l) 

DO 506 1=1, CT 

Fca#i)=rco,i)-vEcici} 

00 506 Ksl,CPl 
eCI»K+C)aBa,f<+C3-ALFA(;i 
CALL INVpHtCB) -• 

CALL BCMULCO) 

CALL BFMULCaS 
CALL APLMULCJZ-I) 

CALL MBAPLCOZl 
CONTINUE - 


,K) 



JZ=66 

CALL ABJSBOT 
CALL TRlAtJJ 
CALL TR10CO) 

CALL APVMULtJ-l) 

CALL AaMULCO-l) 

DO 507 1=1, -CT 
F(a,I)=FCJ,I)-VEClCl) 

DO 507 K=1,CP1 
B(I,R+C)=BCl,K+C)-ALFAti,K) 

DO 511 1=1 ;CT 
aETAiCI3=0(10,I) 

DO 511 K=1,CP1 
aETA2CI»K}sP(56,I#KJ 
DO SOB KKK=9,l,-i 
K=KKK 

CALL PVPLCK) 

DO 509 1=1, CT 
DO 509 L=1.CP1 

6ETA2CI;Ll=FClZl+K-l,IrL)-ALFA(l,L) 
CALL PVUNOL 
DO 510 1=1, CT 
aETAlCIJ=OlK,I)«VEClCI) 

CONTINUE - 


CALL 

CALL CPLMOHU COFF , »ETA2 , ALFA) 
CALL CQMULCO) 

DO 512 1=1 ;CT 

S?I,K+C)=B|l£K+C)-ALFAtl,K) 
CALL INVPK’hB) ^ 

CALL APLHUHaZ-l) 


CALL BFMULCJ) 



542 




544 

543 

I mmmnm: 


545 

546 

548 

547 

54f 


551 

550 


CALL AQ«ULCa«l) 

DO 542 1=1, -CT 

f CU,I)sFcg,I)-VEClCI) 

DO 542 K=1 ,CP1 
Ba»K+C)=BCJ K+C3»ALFAC1,K3 
CALL I/JVPHTIB) 

CALL BFMULCJ) 

CALL APL^«U}ilDZ-l3 
CALL MBApLCaZ) 


J=3b 

CALL CONBCCJ.O) 

CALL INVPHTCB) 

CALL BCMULCJ) 

CALL BFMULCil) 

DO 543 ua=37,39 

jRja 

CALL TRIACJ) 

CALL TRlBtJl 
CALL TRICCJ) 

CALL APV«0L(J«1) 

CALL AQ«UL(k3-l) 

00 544 1=1, CT 

Fca,i)sF(a,i)-vBCiCi3 

00 544 K=l,CPl 
B(I,K+C)sBa,K+CJ*ALFA(l,K) 
CALL BCNULfuS 
CALL BFMyLCO) 
eONtlNUE ' 

J=40 

az=az+i 

CALL TRIAta) 

CALL TRIBCJl 
CALL TRlCtal 
CALL APVMULU-i) 

CALL AaMUL(i3»l) 

00 545 1=1 ;CT 
FCa,l)=FCa,l)-VEClC13 
DO 545 K=1,CP1 
8ClfR+C)=BCI,K4C)-ALFA(l,K3 
CALL INVPK'ft&J 
CALL BCMULCU) 

DO 546 1=1 ;CT 
6El'AlCI3=al35#I> 

00 547 KKKa34, 21,-1 
KSKRK 

CALL PV«?1UL 

DO 548 1=1, CT 

8EXAltIJ=«(K,I)-VKCltl} 

CON'lINUe - 

DO 549 1=1, CT 

DO 549 K=1,CP1 

BEtrA2CI,KJ=FC79flfKl 

DO 550 KKKc34, 23,-1 

K=KKK 

KZ=79-K 

CALL PVFLCK3 

00 551 1=1, CT 

iin S*il .List €Pi ■ 

aETA2CI,LJ=P(KZ,I,LJ-ALFACl,L) 

CONTINUE 

CALL P¥PV(2l,223 

CALL OFFCtJ,21,Tt213) 

CALL CPLMOLUrCOFF,BETA2,ALFA> 
CALL COBULCD) 

DO 552 1=1, CT 



552 

553 

554 
4 

r 

# 

555 

» mm< 


557 


555 

I mm 


558 

4 

I 

«a® IP 


rCj,I)5FCU,I3-VEClC43 
CALL. MBAPLCJZ) ^ 
CALL. BCMUL(U) 

CALL BFMOLCU) 


3=41 

JZS81 

CALL TRIALJ) 

CALL TRJtBtJ) 

CALL TRJtCtJ) 

CALL 0rrALa,20.T(20Ji 
CALL ApyMOLta-l) 

CALL AQMULCa-13 

DO 553 1=1, Cl 

FCJ. Dared, i)-VEClCl3 

DO 553 K=1,CP1 

aCI,K+C)aBCl,K+C)-ALFAtA,K3 

CALL INVPKTCB) * 

CALL BCMULCd) 

CALL AQMUL(d-D 

DO 554 1=1 ;CT 

FCd, Dared, I)-VEC1CD 

CALL BFMULtd) 

rouitipiy cc=A.P{8y) 

rouxtipiy alfa = Aorr.pcbe) 

DO 555 1=1, CT - 

DO 555 K=1,CP1 

« h f I ' + ALF A C I , K ) 

CALL'MBAPLCdZ) 

‘ IIP «4k » «itt Ml M IM III «p fW W w «K '>«l «!» in» 

JZ=81 

DO 556 dja42,45 
d»da* " ■ 

dZsJZFi 
CALL fRlklO) 

CALL TRlBLd) 

CALL TRICtJ) 

CALL APVMBLtd-D 

CALL AQMULCd-D 

DO 557 1=1, -CT 

F(d/DarCd,D«VEClCD 

DO 557 K=1,CP1 - 

B(I,K-»-C)=tt(I,K+C)-ALFAtl,KJ 

CALL INVPRTCI&) 

CALL BCMULCd) 

CALL BFMULCd) 

CALL APLMULCdZ-l) 

CALL MSAPLCdZ) 

COWTINUe • 


CALL 

CALL 

CALL 

CALL 


d=46 

dZsHS 

CALL TRlAtd) 

CALL TRiBCd) 

CALL TRietd) 

CALL APVMOLIJ-D 

CALL AOMULCd-D 

DO 558 1=1, CT 

FCd/DsrCd;i)-VEClCD 

DO 558 Kai.CPl ^ * 

B(I,K'f*C)=BU,K+C)-ALFACl,K) 

CALL INVPKlCfi) 

CALL BFMULCd) 

CALL APLMULtdZ-D 
subtraction CCaCC-ALFA 
CALL mtMULZCB.CC) 
multiply Pt46)«B#CC 


12 






ds47 

CAijL REBABCd.O) 

CALL APVMUHJ-l) 

CALL AQMUL(OI-l) 

DO 564 l=l;CT 
FCj,I)crC‘J,I)-VEClCI) 

DO S64 Kisl «Cpi 

t >64 8 Cl/K+CJ=tttl,K 4 C)**ALFA(l,K 3 

CALL BFMULCJ) ' 

8ACK«SUB6TJlfUtiOBi 

HB ««•••« m M M ■*•*«••• W K ■■• w W ■• M w •* aa « w a 

a=46 

CALL OHPVXCa) 
dZsHS ■ 

DO 570 da=45,40,*l, 
jaja- 

CALL OMPVXCJ) 

CALL QHPLXCa,dZr473 
JZ=dZ-l - 

570 CONTINUE 

00 571 UJ» 39 , 36 r»l 
jsdd- 

CALL OHPVXCd) 

571 CONTINUE ' 

dZ=79 

CALL Q«PLX(a,JZ , 473 
DO 572 aiJ=34j23,*l 
JsaU' 

(JZbJZ -1 

CALL 0MPVXfJ3 

CALL QMPLX <0^02,473 

572 CONTINUE ' 

JZa67 

DO 573 Jd* 22 , 21 ,-i 
jsao' 

CALL GtMPVXCiJ3 

573 CONTINUE 

az=66 

CALL QMFLXta,JZ#473 

DO 574 

a=jd 

JZ=JZ»l 

CALL QHFVXCd) 

CALL QMPLXC«J!,aZ,473 

574 CONTINUE ' 

171*20 

JZ*56 

CALL QMPLXCd,aZ,203 

DO 575 aj=9,i,-l 

JsJJ 

JZ*JZ-1 

CALL OHPVXCd) 

CALL QMPLXCvI,az,203 

575 CONTINUE ' 


STOP 

END 

SUBROUTINE ENVCTT»HV3 

dimension HV€Z 3 

riV(13=l. 

riVt23*l. 

RETURN 

END 

SUBROUTINE ENLCTTrHL3 
dimension HL (23 
HLC1)*1. 

HLC23=l* 




KETUKN 
£ND 

SlIBROtlTlNE riNDKCTX,AK) 

OIMENSIUn AKt2) 

AK(.U = 1, 

AKC2)=i. 

RETURN 
END 

^ Wft <!s® fp *P9 sst pBi HP IP IP m IP P* IP PI PS PI wp ®i pi ifp ppi HP nii pi p ip p -p pp up ig^ pui pp tip pp ip' ip fp up m; pp pp pp ^ ,pp pp pp .p^ pp p^ ^ 

SUBRUUfANE CONBCCU, J.TYPE) 

•SUBROUTINE fOR CUMRUTING THE ELEMENTS Or "B" & "C* SUBRATRICES FOR 
tA PAR'i’fAL MH TOTAL FOR PARTIAL , mPE^I FUR TOfAL 

INIEGLR L f C“i » J, wOL , CT * - 

COMMON N.C,CPi,TKi|lJC/CT 
COMMON /XB/BC9.y)/XC/CCt9,9) 

COMMON /XVAK/5Lt47,4J,BIiC47),RLC47) ,SV(47 ,4) , BV (47 ) ,RV (47 3 , T C47 ) 
COMMON /XCONST/ETAC473,»SLC47),SSVt47} ' V ' . . 

DIMENSION XRt4,43 ,TR(4,43,AK£4) ' - 

DO S Ia|,CT - 
DO 5 Ksl.CT 
5 CC£1,K3=5. 

GLS-C1.TSBLCJ33 

GV»-Cl.+SSVCa)3 

E8XLaETA(J)/BLCJ) 

DO ID 1=1,0' 

CC(I,I+C)aRV(J+l3 

8CI»CT)=0» 

B(I,1)sGL 

10 BCl.i+OsGV 
IF£IT5fPE*EQ,l)GDTO 100 
CALL riNDR(Tca3#AR) 

CALL DKBXDX(T£a5,XK3 
CALL OKBXD¥(TCtI3,llO 

DO 20 1=1 »C 

El=-SL(D,I3=AK(I3/BLCa3 
£2»SVC«3,r3/BVtO) * 

DO 20 K=1,C 

BCI+C,K3=E1 

8CI+C,K+C3=B2 

irtliME.KlGOTO 11 

B(1+G,K)=BCI+C,K3+AKCI3 

BCI+C,K+C3aBCI+C,K+C3*l. 

11 SUMX=a, 

simifso, 

DO 30 IFsl,C 
DELFK=0, 

IF(K.Ea.IP)DELPKal, 

SU«X=SUMX+XR£I,XP3*tDELPK-SLCa,IP>/BL(a33 
30 SUM¥=SUM¥TXKCI,IP3*CDELPK-S¥Ca,iP3/BV(D)3 

BCl+C,K)=EB¥L=l£SLt J,-|3»SUMX/aL£a)4-8(3:+C,K3 3 
20 SCI+CfK+OalEBfL^SLiJiIilSUMr+iClTlC^K+Cn/fiVCD) 

CALL DKB¥UT(TCa),AKJ ' ' . 

00 70 1=1, C 

70 8CI+C,Cf3=E8YLI5Lia,l3IAK£I3 

ir(ETACJ),EU.l,)RKTURW 
FAClsCS,-ETA(J3)/tBVtJ+i)=lB¥CJ-l-l3 3 
DO 80 1=1, C • ■ 

DO 80 K=l,C 
DELIK=0, 

IFCI,EQ.KJDELIK=1, , » 

80 CCCItC,KTC3=FACl»(DELIK*^BV£a+13*SVCJ4l,I)) 

RETURN 


t 1/m mmmn 




IFOR TOTAL CONDENSERS 
too S0V=l./lBVCO3IBV(a)3 

S0L=l./tBLCU)=BL£D)3 
DO 105 1=1 ;C 
El=-SLCU,l5*SaL 



102 

105 


E2=SVCJ,IJfSQV 

DO 105 KbJ,-C 

IF( I,EQ,K1GOTO 102 

BCl+C,MsEl 

ti(I+C,K+C)=:E2 

GOTO I 05 

BCI+C,KJ = tBOCa)-Sli(d,I)|»SQL 
BCI+C,K+CJ=CSVca,X3-BVCa))*SQV 
COHTIfiUE ' • 

RETURN 
END ■ 


SUBROUTINE KEBABCU - ITXPE) 
SUBROUTINE FOR CUMFUTING THE ELEMENTS OF 


» IP « I)* in iw m MF «• 


_ & *^8**’ SUB WAf RICES FOR- 

lA PARTIAL OR TOTAL REBOILEKJIUPESO FOR PARTIAL , ITTPE=1 FOR TOTAL 
INTEGER C#CPl,TNDC,CT 


10 


20 


2A 


27 


30 


40 

50 


V,# V,)ri f iWUV.VA 

COMMON N,C,CP1,T«0C,CT 
COMMON /XA/AC9!9)/XB/B(« 




V n Sufi'S f / 1..^ V ^ i 

COMMON /XVAR/SLC47,41,BLt473,RLt47),SV(47p4),8V(47),RVC47),Tt47) 
COMMON /XC0NST/ETAf471,SSL(47),SSV{47) * 

DIMENSION AKC4),XAC4,4)#yKC4,45 ' ' 

GVa^d^+SSVtvJ)) 

GL=-C1,+SSLU)) 

DO 10 1=1 »CT 
DO 10 K=lrCT 
ACItRlaO, 

DO 20 1=1, C 
A(I,X)=RLCJ-1) 

B(I,I)*GL- 

8CI,I+C3=GV 

BCI,CT)=0. 

IFtITXPE,EQ,l)GOTU 50 
CALL FINDKCTCJ3,AK3 
CALL DKB¥DXtT(a),XK) 

CALL DKByDytXCa3,5(KJ 
DO 30 1=1, C 

E1=»SLCJ.1)^AK(I3/BLCJ) 

E2*SVCd,I3/BVCd) ■ 

DO 30 K=1,C 

B(I+C,K3=E1 

B(i+C-K+C3=E2 

IFCI,NE.K3G0T0 24 

BCI+C,K3=BCI’»'C,K}+ARC13 

B(I-!'C,K+C)=BCI+C,K+CJ"1. 

SUMX=0, 

SUMY=0, 

DO 21 IF=1,C 

DELPKaO. 

IFCIP-E«.KJDELPK=1, 

SUMX=SUMX+lRCI,XPJ#tDELPK-SL(J,IP)/BLCJ)) 

SOM¥»SUMY+yKCI,XP)=tOELPK-SVCd,IP3/BVCJ)5, 

B(I+C,K) = tSLCa,I)»SUMX/BLCJJ+BCI+C,Kn/BLCa3 ,, 

8CItC^K+C3slSLtd,l)»SUMl/BLCd3+BCI+C,K+C3J/BV(a) 

CONTINUE ' 

CALL DKBYUriTCa3,AK3 
DO 40 1=1, C ^ ^ 

B(I+C,CT)=5LCa,I)=AK(I3^BLCJ) 

RETURN 

SQL=l./tBLCU3»BLCd)3 

SQV»l./tBVCU)*BVCU33 
00 60 1=1,0 
El=-SL(U,13»SaL 
E2=SV(J,I3*SQV 
DO 60 K=i,C 
IFCI.EQ,K3GUTa 51 
B(1+C,K3=E1 
BCH-C,K+C3sE2 
GOTO bO 



1ft 


bl BCI + C,KJ=!CBIiCU 3-SIL«CJ,1:) IfSQL 

tiCi+c,K+c;=tsvta,i3-§vtU))»SQV 
60 CONTINUE . ' " v / / w 

RETURN 
END 

subroutine TRIAIUI 

iSUBROUTiNE TO COMPUTE THE ELEMENTS OF THE TRIDIAGONAL «A" MATRICES 
INTEGER CrCFl .TWOC-CT - ' 

COMMON N,C,CP1,TWUC,CT 

COMMON /XA/AC9-y)/XC0N5T/ETAC47) ,SSL(47),SSV(47) 

dimension HLC4)rDHLl43 ' ' - - 

DO 10 1=1, CT -- 

DO 10 K=1,CT 
ID ACI.KIaO, 

CALL ENLCTta-i; 

CALL ENVCrC» 

SlfM=0. - - 
DO 12 1=1, C 
ACI,13sHLU-l) 

ACCT,I)=RLCU-1)»HLC13 
12 SUM=SUM+SLC0-l,I)4DHLtIJ 

ACCT,Crj=K£iC0-13»SUM - 
RETURN 


,a-13,HL) 

!u«U,dml) 


I' 


END 


timmm:mmmmmmmmmmmmm<mmmmmmmmm» '. 


SUBROUTINE TRIBCJl 
ISUBROUTINE TO COMPUTE ELEMENTS OF THE TRIDIAGQNAL ’’B" SUBMATRICES 
INTEGER C,CF1,TW0C.CT 

COMMON /XB/HC9-9)/XC0NST/ETAC47),SSLC47),S5VC473 
COMMON /XVAR/Stt47,4),8L€47),-RLC47) ,SVC47,4)/BV(473,RVC47),TC473 
DIMENSION AK(43,XRC4,4),TK(4,43,'HLC43,HVC4>,0HLC4),0MV(4) 
Sl=-U,+SSLCa)3 " ' 

S2s-Cl,+5SVCJ)3 
DO 10 1=1, C 
8CI,I3=Sl 
8CI,i'»"C3sS2 
10 8(I,CT)=0, 

EBTL=ETA(U3/BL(U3 

C0EFF=E8TL/BL(J) 

CALL FINDRCTCU3,AR3 
CALL OKBXDXlTCd>,XK) 

CALL DKSX0YtTCa3,IK) 

DO 20 1=1, C 
E1=-AKCI3’^5L(J,I3 
S2»SVCU,I3/BVCa3 
DO 20 K=1,C 
0(I+C,R3=E1 
B(ITC,K+C3=E2 
IFiI,NE.K360T0 21 
B(1+C,K3=BCI+C,K3+AK(I3»8L(J3 
B(I4e,K+C3sB(I+C,K+C3-l, 


21 

iir 

I ■' 


24 

20 


30 


DV 4 .TV>,t 

SUMXaO, 

5UMT=0, 

’•KC13" VALUES ARE COMPPSITIUN •INDEPENDENT THEN REMOVE THE "I 

• - GOTO 24 - 

DO 24 IP=1,C 
DELPK=0, 

l5SSSji8sJ«i;pr![pUcD.L.K-|«a,i|,jBL(an 

SOMy=SUMY+IKaJlPJMtDELPK-SVCa,iP348VCa33 
CONTINUE . . , 

a(l+C,K3sCO£FF'»:CSL(wJ,I3*SUMXl.8(I+C,K33, , 

B(I+clK+C3 = lEBIL4SLia,lJ»SUM5f4Bt;I+C,K4C33/BV(J3 
CALL DKBYDTtTC J3 ,RK3 
DO 30 1=1, C 

BCl+C,CT3=EBYL45Liafl3»AKCl3 



17 


33 




CALIi ENliCTCJ) ,H1U) 

CALI^ ENVC1’(U),Hv5 
CAtili DHUDXCTCUJfDHti) 

CAl4lj OrtVDT(T(J),D«V) 
S(IMX=0, 

SUMlsO. 

DO 33 1=1. C 

a(CT,I3=Sl»H|4(I3 

8(:CT,I+C)=S2*HVCIJ 

SUKX=SU«X+SL<ca,I3*DHtiCI3 

SUMY»SU«y+5V(J,l)*DHV£|3 

8CCT,CT3=Sl*SOMX+S2»SOMy 

REfURN 

END ■ 


SlIBROUflNE TR1C(J3 
iSUBROOTINE TTO CUMfUTE tHE 


i III fir « if* , 

ELEMENTS OF THE TRIDIAGONAI, "C* MATRICES 


10 


11 


12 


15 


15 

20 


22 

21 

30 


32 

31 


INTEGER CfCPl.TWOC-CT 
COMMON N-C,Cpi,TlfiOC,CT 

COMMON /XC/CCC§,9)/XCONST/ETACf7)#SSL(47),SSVC473 

COMMON /XVAR/S0(47-4),Bt»£47),Rtit473£SV£47,43,BVC47),RV(47),TC473 
COMMON /Xi.INK/0ENOftl473#N0FF,IX(63, 11(63, t.ORVC&3,RATiOC53 
DIMENSION VECT0R(4)^HVC4),DHV£43 ' 

DO 10 1=1 #CT ^ ■ 

DO 10 K=l,CT 
CCCI,K3*0, 

CAUD ENV(TCv3+13»HV3 
CAOO DHVDT(T(J-I-13»DHV3 
SUM=0, 

DO 11 I=i#C 

CCCIj,ItC3=RVCd+13 

CCCCT,I4C3=KVCJ+13*HV(I3 

SU«=SUM-t>S¥ (d+1 , 1 3 *DH V C I J 

CCCCf,CX3=RVCa+13*»'SOM 

IF C A8S C 1 , -eta ( U 3 3 . LE ,1 ,6-6 ) RETURN 

DNK=0 

DO 12 1=1 »C 

VECTORCI)=RV(d+13*SV(d+l,I3 
DO 20 KRsl ,WpFr ■ 

IFCCd.EW.xicCRK) ),ANO, tUUR¥CKK3.EO,l3 3GOTO 15 

GOTO 20 

UNK=1 

0L=IYCKK3 

P 0 %-i^ ' " X ® C 

¥ECT0RCI3 = VECT0R(13+RAT10CDL3»S¥(Dl4,I3 
CONTINUE ■ ■ 

IFtUNK.EQ.UGOTO 30 

COEFF=il.-ETACd33/(BV(d+13»B¥Cd4l33 
DO 21 1=1, C ' ■ 

£0=-C0EFF»SVed+l.l3 

£l=COErF4(BV(d+l3-SV(J+l,I3> 

DO 21 R=i;c 
IFa,EQ.K3G0T0 22 
CCaTC,K+C3=E0 
GOTO 21 

CC<I+C,K+C3=E1 

CONTINUE 

RE*rOR H 

COEfF=(l,-ETAC J33*RVCO+l)/CDENOMCd3»OENOM£d33 

DO 31 1=1 ,C , . . 

E0=-C0EKF^VECT0H(13 
ElsCOErF»(D6NOMCd3-VECTUR(I) 3 
DO 31 K=1^C 
IFCI.E0.K3G0T0 32 
CCCI-i*C-K+C)=EO 
GOTO- 31 

CC(I+C,K+C3=E1 

CONTINUE 
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I mmmrn 


RETURN 

END 


SUBROUTINE OFfACU rJEXIT^TT) 

i?!!A5S^S'llR5?i!igE‘’[ir.JSS:sSkHK^I!Ei? wf-diacdnal -a- matrices 

t *** only for liiauid-interilnKS 


J 

UEXIT 

TT 


ff C473 


10 


20 

800 

30 

801 


40 


50 


70 


J IPt 


Stage no, into wnicn tne interiinK-streaffl is entering 
. stage no, from wnicn tne interiinK-stream is leaving 
: temperature of tne interilnK-stream entering stage-j 
CTT may be equal to that of tne leaving stage-uEXlT 
•or IT may be specified) . = 

INTEGER C#CF1,TW0C-CT 
COMMON N,C,CPl,T«OCfCT 
COMMON /XA0FF/A0rrC9,9) 

COMMON /XblNK/DEftiUMt47),N0FF,IXt6)Al¥Cfe),I»0RVC6),RAT10(6) 
COMMON /XVAK/50C4V,4J,Bb(47),RbC475#SVt47,43vBV(47),KVC473 
DIMENSION Ht.C43,DMU43 - V - 

00 iO 1=1, CT 
DO 10 K®l,Cf 
A0FFtI,KJ=0, 

00 20 KK=l-N0rF 

IFC(a,EU,lXtKR) J,AND.CaEXIT,EQ,I¥CKK)))GOTO 30 
CONTINUE 

TYPE 900-JEXIT,a 

FORMATC IX, 'There is no stream from stage', 12,' to stage', 12) 
STOP 

IF(L0RV(KX3,EQ.0)G0T0 40 
TYPE 901# U 

F08MAT(ix, 'subroutine offa is strictly for liquid interlines 

1 and should not be called for vapor streams. Check stage', 123 

STOP - . 

CALO ENOCTTfHli) 

OO 50 1=1, C 

A0FFCI,i)=RATl0tKK3 

AOFFCCT,l3»RATIO(FK3*Ht,(l3 

IFCABSCTT-TUEXIT3).GT,1.E*63 RETURN 

GAbti DMLDT(TT,DHL.3 

DO 70 1=1, C 

AOFF C CT , CT 3 = AOFF C CT , CT 3 +Sli ( JBX IT , 1 ) »DHb C 1 3 

ADFFCCT,CT3=A0FFCCT,CT)»RATI0CKK3 

RETURN 


END 


» «l* fW P» IM M « P» IW W lli «■» Pi «m> » «» W N 


» «, l» MM «W w «» ■» fW « 


*' SUBRUUTXNiii Ot’^FCtU kJI^XXT TT'l 

iSUBROUTINE TO COMPUTE ThI: EUEmIiNTS OF THE OFF-DIAGONAD "C* MATRICES 
ICHAt/lNG structure UF "C-^SUBMATRICES) 
f*** Only for Vapor-InterilnKs *** 


: JEXIT 
J Tf ■ 


10 


20 

900 


Stage no, into which tne interllnlc-streain is entering 
j stage no, from wnicn tne interilnx-stream is leaving 
s temperature of tne interiinX-^stream entering stage»a 
ITT may be equal to tnat of tne leaving stage»JEXIT 
*or TT may be specified) - * 

INTEGER C,CPl,TWOC-CI 
COMMON N,C‘ Cpi,T«0C,CT 

COMMON /XVAR/SbC47!-45!BU(47 3,RLC47>,5VC47,43,BVC47 3,R¥C47 3,T(47) 

COMMON /Xl,INK/DENOAc4n,bFFjlX<6),iYCe3,f,ORn63,RWlO(63 

dimension HV(4),DHyC43,VECTOR{43 

DO 10 1=1, CT 

DO 10 K=1,CT 

C0FFCI.F)=0, 

S n li K K *** 1 w f 1 F* 

IFCC J»£u7li( IKK) 3 .AND, CaEXIT.EO.lYCKK) 3)60TO 30 
CONTINUE ’ 

TYPE 900,UEXIT,U ^ to , 

FORMATC IX, 'There is no stream from stage', 12, 

STOP 


to stage.- ',123 
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IF(I<0RV1KK),E0,1)WQ 4U 
TYPE 901, J 

FORMATUX, 'Subroutine OFFC Is strlctlY for vaoor interlinics 

sbouid not foe called for liquid streams. Check: stage', i2) 

STOP 

CAI.L EflV(TI,HVl 

DO 50 1=1, C 

COFFCI,I+C)=RATIOtKKl 

C0rFCCt,I+Cl=RAfI0CKK>»HV(I3 

lFCABSH'T-tia£XIIl),Gr,l,E-S3G0T0 70 

CAE.Ii DRVDf(Tf,DHy3 

DO 60 IBI.C 

COFFtCT,CXJ=caFrtCT,CT>+SVCJEXlf,l3skDHVCl3 
C0FFCCT,CX)=RAT10tKK3»CUFF(CT,CX) 
IFCABS(l,-E'rA{a3).tife:.l.E-6JRETURN 
DO -80 1=1, C 

VECT0RCl)=RV(j4l)*SV(vJ+l,I) 

DO 90 Dfo=l,NaFr 

IFCta,E«.IXlDD)).AND.CLORVCLD).EQ.J))GOfO 95 
GOTO 90 ■ ■■ ■ ■ ■ 

DO 99 1=1, C 

yECT0Ra) = VECXURCl)+RAT10(H.)»SYCIY(LL),I) 

CONTINUE 

FAClsCl.»ETACa33*RAT10(XK3/DENa«(JJ 
DO iOO 1=1, C ■ ' 

EOa-FACl»VECTORCIl/DEND«Ca) 

El*rACl*(l,-VECXUKC13/DEN0M(Jn 
DO 100 K=1,C 
IFCI.eQ*KJGOTO 110 

corrci+c,Kf€3=Eo 
GOTO *100 

corra+c-K+C)=Ei 

continue' 

RETURN 

END 


iXhls sworoucine compures me t«sc row or one 'w- BUDwat* 

'depending on tne type of speciticatlon variafoXecDABEli) for any stage a 

j I stage - * . 

ITYPE : if* stage U is a condenser or a reboiler 

then set itype = o tor a "partial" condenser or reboiler 
else assign any integer to iTf PE ^ 

YAD s if no specif ication-vaiue is required 

then assign any real number^ to VAD ^ ^ 

else assign given value of the specification variable to VAL 
ICOM : component no. tit not required then assign any integer! 

INTEGER C,CPi,TWDC,CI 
COMMON N,C,Cpi,TauC,CT 
CDMNON /XA/A(9,9)/X8/BC9.9)/XC/CCC9,93 

COMMON /XVAK/s{.C47,4),Bt.C47),RLC47),SVC47,43,BVt473,RVC47 3 ,Tt47 3 
COMMON /XC0NST/ETAt475.SSGC4?),SSyt47);^^,,,^- ' ' 

DIMENSION HV(4),Hyplt4),HLC43,fiDMU43|DHy(43,DHD(4) 

GOTOtlO,2«, lO,4U,b0^6O, /0,80,90,IUO, 110,120,1303LABEL 

TYPe-90a, J,I iABED ' ^ ^ Til 

30 FORMATCix, terror m specification, for stage',13, 'spec, type ,133 

STOP" ‘ '* ' " ’ ” 

NO CONDENSER SPECIFICATIONS 
3 Sl»»tl,*SSDtiJ3) 

S2=-Cl,-t-SSVtJ3) _ . 

iFCIfYP£,EO.l)GOTU 15 
CAhL EHYtTCOl.HV! _ 

CADD ENYtTCO+l) ,HVpl) 

CALL DHVDTCTCU3 ,DMV3 
SUM=0, 

DO 12 1=1, C 
BCCT,I3=S1*HVPHI3 
BCCT,I’fC)=S2*HVCI3 


SUBROUTINE PROBtU, LABEL, ITYPE, VAL, ICOM) 
subroutine computes the last row of the "A", 


submatrices, 
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2 SUK=SUM+SV(U,I)»DHV(1) 

B(CT,CTJ=S2*SUM 

CALL OHVDTCTCJ+UfDMV) 
cccca’,cT)=o, 

DD-13 1=1, C 
CCCCT,I)=0, 

CCCCT,I+CJ=KV(J+lJi«'HVPm) 

3 CCCCT^CX)=CCCCT,CT) + (KVU+l3fSVta+l,I)+Sj»SLW,l))»DB¥(I3 

,5 CALL ENLCT(a3,HL) 

CALL 0HLDTCT(J3,LHL) 

B(CT,CT3=y, 

00 16 1=1, C 
BCCT,I)=S1»HL(I) 
afcr,i+c)=s2»HLCiJ 

8CCT,CT3=:B(Cf,CT) + CKV(J+l)»Sy (0+1,1 3+S2»SVCJ, 13+51 fSLia, I) 

1 0HLCX3 
ec(Cf,i3=o, 

.6 CCCCT,I+C)=KVC0+U+Kyia+1)+HLCI) 

CCCCT,CT3=0, 

REfORN 

NO HEBOILER SPECIFICATIONS 
!0 Sl=-Cl,+SSL(a33 

S2**Cl.+SSV(a33 
IFCITTPE,EQii3G0T0 25 
CALL EML(T(a-13,HLMl3 


CALL EfiLCT(a3,HL) 

CALL DHLDXCTCJ-13 ,0HL3 
ACCT,CT3=0, - 

00 22 I»lrC 

ACCT,I3orC»CO-13*HLMI (13 

A(CTjCT3sA(CTjCT3+CKLta-l)+SL(0-l,l)+S2+SV(a,I33*DHLCI3 
BCCT,I3=Sl+HLa3 ■ 

22 &(CT,I+C)=S2+HL«1 CI3 
CALL 0{4LDT(T(J3,D«L) 

SUH=0, 

DO 23 1 = 1, C 

23 SUM=SUM+SL(U,I3»DHLtI3 
B(CT,CT3=St+S0M 
RETURN 

25 CALL BNL(T(d-13,HLMl) 

CALL DHLDTCTCJ-13, DHL) 

A(CT,CT)=U. 

DO 2b 1=1, C 

ACCT-I3=KL(U-.13+HLR1CI3 

A ( CT I i 

A(CT;CT)=A(CT,CT) + (RHa-13 + SLCa*l,I)+Sl»SLCa,I3+S2»SVCJ,I)3* 

1 0HL(I3 - \ ■ - 

B(CT,I3=SI+HLM1113 

26 B(CT,I+C3=S2»HLM11I> 
aCCT,CT3=0, 

RETURN 

« CONDENSER-HEAT-DUTy SPECIFIED 

30 Sl=-Cl*+5SL(a)3 
S2=-C1*+55V(J33 
CALL'ENytf(U+l3,HV3 
CALL DHVDT(TCUtl3,DHy3 
SUM=0. 

DO 31 1=1, C 
CC(CT,I3=0, 

ccccT,i+c3sKvca+i)+Hy(ii 

31 SUK=SUM+5VCd+l,I3+DHytI3 
CC(CT,CT3=RV(J+1)*SUM 


DHL3 


♦DHytll 

*SUM 


V V VV .Jl I V * ^ \ ^ W rj 

IF(ITyPE,EO,l)GOTO 35 
CALL EflL(Ttd3,HL3 
CALL ENV(TCJ3,HV3 
CALL DHLDTCT<a3 ,DHLJ 
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35 


36 


CALL DHVDtlf C^JrOHVJ 

£1*0 (I 
£2*0* 

00 32 I=lrC 

SCCT,I)=51*HL(I) 

tJCCT,I+C)sS2»HVa) 

ei*Ei+sLca,i;)*DHLti) 

E2=E2+SVC0-3:)»DHVCI) 

B(CT^CT3*6i»EJ+S2»E2 

CALL ENLCy(iJ3,HL) 
CALL DHLDtC‘i;Ca),0HL) 
£1=0, ■ 

£2*0*. 

DO 36 1*1, C 
BfCT,I)*S14:ML(I3 
B(eT,I + C)*S2>»=HLCI3 
ei*Ei+SLC0,l)*DHLtl3 


E2*E2-l-SVCJ,JC)*0HLtI3 
BCC'r/CT)=Si*El+S2*E2 
REfURN 

I REBOILER-HEAl-DUTY SPECIFIED 

40 Sl=-Cl,+SSLtJ)) 
S2»-C1.+S6VCJ)3 
CALL-ENL(Tta-l5fHLMl) 
CALL DHLDTCTCJ«J)#DHL) 
SUM*0* 

DO 41 1*1, C 

ACCT,I)aRL(J-13*HLMl(I) 

ACCr,H-C3=0. 

41 SU««SUM+5HD-1,I)»DHLCI3 
ACCT,CT3*KL(J-134SyM 
IFtlWE.EQilJGOia 45 
CALL EMLC‘r(a),HL) 

CALL eNVCT(a),HV) 

CALL OHLOf (T(J),6 hL) 

CALL Di4VDXCfCD)»D«V) 
El*0, 

£2*0, 

DO 42 1*1, C 
B(Cf,I)=Sl»HLCI) 
B{CT,H'C3=S2=*'HVtl3 
El*El+SLCa,I3*DHLtI) 

42 E2*E2+SV(J.I)*DtlVCI) 
BCCf ,CT)a5i»El+S2»E2 
RETURN 

45 CALL ENV(TCU),HV) 

CALL 0HVDTCTCa),DHV3 
SUM=0, 

DO 46 1*1,0 

S(CT,I3aS14HLNlCI3 

0{CTIi+C)=S2*HVCI> 

46 SUMsSUM+SV(J-I)*DHVCI) 
8CCT,CT)*S2*SUM 
RETURN 

i REFLUX RATIO IL/DI SPECIFIED 

50 DO St 1*1, C 
a(CT,I3*l, 

51 BCCT,I+C)e-VAL 

gSls'fii £t 

52 CCiCT,I)*0, 

RETURN 

i REBOILER RATIO SPECIFIED 

60 DO 61 1=1, C. 
6CCT,I)*-VAL 

61 B(CT,I+C)=1, 

rTi=f» 
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i T(3) 
70 


i vcu 
80 


i 111 C J J 
fO 


specified for re&oiiers) 


I sv (3 
too - 


i Sblj 
XIV * 


i yUf 
120 ' 


dCi) specified for condensers? 


1 x( 3 , 

130 - 


DO 62 1=1, CT 
ACCT,I)=0, 

RETURN 

SPECIFIED 
DO 71 1=1, CT 
ACCT,I)=0, 

B(CT,I)=0. 

CCCCt.I3=0. 

dCCX.CT3=l, 

return 

SPECIFIED Cor D specified for condensers) 

00 81 Is1,CT 
A(CT,I)=!0, 

CCtCT,I)=0, 

00-82 1=1, C 

0CCT,I)=O, 

accT,i+C 3 =i. 

B(CT/CT)=0, 

RETURN 

SPECIFIED lor B Specified for reooilers) 

DO 91 1=1, CT 
A(CT,I)=0, 

CCCCT,I)=D, 

DO-92 1=1, C 
BCCT,I)=l, 

B(CT,I+C)=0* 

BCCT,CT)=y, 

RETURN 

,1) SPECIFIED (or dli) Specified for condensers) 
DO 101 1=1, CT - 
ACCT,I3=0, 

BCCT,I)=0, 

CCICT,I)=0, 

BCCT,C+1C0H)=1, 

RETURN 

,1) SPECIFIED cor blij Specified) 

DO 111 1=1, CT 
ACCT,I)=0/ 

BCCT£l)=D. 

CCICT,I)=0. 

B(CT,iCOM)=l, 

^ m l| O fej 

1) SPECIFIED Cor xdU) specified for condensers: 
DO 121 1=1, CT 
ACCT.IpD, 

CCCCT,1)=0. 

D0'122 1=1, C 
B(CT,I)=0, 

B(CT,I+C)=»VAD 
a(CTlc-HCO«)s;BCCT,C+IC0«3 + l, 

B{CT,CT3=0, 

RFTUHM 

13 SPECIFIED lor xbli) specified for rebollers) 
DO 111 1=1, CT 
ACCTiI)=0. 

CClCf,I)=0, 

00132 l=i,C 
6{CTiI)=-VAD 
BCCT,H-C)=0, 

6CCT,IC0M)=BCCT,IC0«)tl. 

END 


bli) specified) 


xdll) specified for condensers) 


subroutine PVPLCaPV) 

COMHQN N,C,CPi,TWUC,CT 

COMMON /iP^PC8§,9,S)/XBETA2/BETA2t9,5)/XADFA/ADFA(9,5) 
DO 10 I=t,CT ■ ■ ■ ■ - ' ■■■' ' 
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UO K=lrCPl 
AI.FA(I,K)=0,0 
L)0 10 L*l,Ct’l 

1 J ^ t ^ J ^ ^ I ' ^ C I^+C , K ) 

HEftIRM ■ ' ■ 

ENO 

SUBROUflNE PVPVCaFVl,aPV23 
integer C,CP1,TM0C,C| 

COMMON N.C,Cpi.TliOCfCf 

COMMON /XP/P(8§,9,55/XAliFA/ALFAC9,53 

DO 10 1=1, c*r 

DO 10 K=l,CPl 

ADFA(I,K)=0,0 

DO 10 l.*l,CPl 

Al,FACI,K)rADFACI,R3+PCJPVl,I,D)»PCaPV2,D+C,K) 

RETURN 

END 

SUBROUTINE PLPV C U1 , 02, Ui) 

INTEGER C,CPl,TWOC,CT " 

COMMON N,C,CP1,TW0C,CT 
DIMENSION Ol(9;5),U2t9,!5),U3C9,5) 

DO 10 1=1, CT ■ * 

DO 10 K=i,CPl 
03a^J|^=UlU,CT)*U2lCT,'K3 

U3|l^K)=oS(I,K)+UI(I,D3»U2a,K> 

END - 

SUBROUTINE PDPDtUl,02,Ui) 

INTEGER C,CP1,TW0C,CT 

COMMON N,C-CPl,TiiOC,CT 

DIMENSION OlC9,5),U2(9,5J,U3(9,§) 

DO 10 1=1, CX- ■■ ^ - ■ 

DO 10 K=l,CPt 
U3CI,vJ)=Ul(I,CT)»U2tCT,A) 

DO'IO 0=1, C 

U3tI,K)=U3CI,K)+UHI,03»U2(D,K3 

RETURN 

END 

SUBROUTINE APVMUDia) 

integer C/CPi,Twoc,cr 

COMMON N,C,epi,TWUC#CT 

COMMON /XA/A(9j9)/XP/P(8S,9,5)/XADFA/ADFAC9,53 

IFCABSCHRR*1.03*DE.1.0E-063 GO TO 30 

DO 10 1=1, C 

DO 10 K»l,CPl 

ADFAC1,K3=RRR»P(J,1,K3 

GO TO 40 

DO 20 1=1, C 

DO 20 K=l,CPl 

ADFAlI,K)=PtJ,I,X3 

DO 50 K=l,CPi 

ADrA(CT,K5=ACCT,CT)»PCa,CT,K3 
DO 50 Li» 1 , C 

ALFACCT,K5=ADFACCT,K3+AtCT,L3»Pta,D,K3 

RETURN 

END 

SUBROUTINE APDMUDU2M13 
INTEGER C,CP1»TW0C,CT 

COMMON N,C,Cpl,TW06,eT ^ 

COMMON /jCA^A(9,9)/XP/PC85,9,53/XADFA/ADFA(9,53 
RRRbA(1«1) " . , 

IfWsC 1.0-RRR).DE,1.E.063 GO TO 20 
00-10 1=1, C 
00 10 K=l,CPl 



24 


10 ALFAtI,K)s=RHR»PCJ2;«i ,I,K) 
GO TO 50 ' ■ 


20 DO 30 I=lrC 

00 30 K=1,CP1 

30 ALFACI.K)s:p(JZMt,l,K 3 

50 DO bO K=lf Ct*l 

A0FA(CT,K5sA(CT»CT?»PCaSSHl,Cr,K) 

DO 60 l.=i,C 

60 ALFA I Cl , K ) sALF AC CT , K 3 tA ICT , 0 ) ’•'PC OZ , 0 , K ) 

RETURN * ' ' ■ ' ' 

END 

subroutine CPVP«Ul,ia/COFF#BETA2,AOFA3 
INTEGER L%CPt,TWOC,CT ' 

COHMON N.C.CPi.TwOCrCT 

DIAENSION CCIFFC9,93,BETA2C9>53,A0rAC9,5) 
RRRSC0FFC1,CP13 

IFCASSClf0-KRR).0E,l,OE-063 GO TO 30 
00-10 I=i,C 
DO 10 K=1»CP1 

10 ALFACI,K3=RHR»BETA2CI+C#K) 

GO TO 40 


DO 20 1=1, C 

DO 20 K=l#CPl 

A£4FACI.K3=BETA2CI+C,KJ 

IF(ABSU,O-ETACa)3*C»E,l,0E-06) GO TO 70 

DO -50 2=CP1,TW0C 

DO 50 K»l,CFl 

AGrACl.K)xo,0 

DO 50 LsCPl/rTWOC 

A£jFACI,K)=ALFACI,F)+COFrCI,G)*BETA2CG,K3 

GO TO §0 - 

DO 60 I=CFl,TWOC 

DO 60 K=lfCPl 

ALFAtI,R)xO,0 

DO 90 K=1»CH1 

ALFACCT,K)xO,0 

DO 90 L=CP1,CT 

A0FACCT,K)=ALFACCT,K)+C0FFCCT,iu)4BETA2(I*,K) 

RETURN 

END 

SUBROUTINE CPBMULiCJ, CUFF, BETA2, AGFA) 

INTEGER CfCPlfTWOC.CT 
COMHON N.C.Cpi ,T«t!C,CT 

DIMENSION COFFC9,y),BErA2C9,S),AliFAC9,53 
RRRsC0FF(l,CP13 

XFCABSC1.0-KRR).t.E,l,i)E-063 GO TO 30 
DO -10 1=1, C 
DO 10 K=1,CP1 

AIirACI,R3=RRR»BETA2(I-^C,K) 

GO TO 40 
00 20 1=1, C 
DO 20 K=1,CP1 
ALFACI,F) = BETA2M + C,K) 

IFCABSU-O-EfACa) J,0£.1,E-06) GO TO 70 
DO -50 I=CP1«TI#0C 
DO 50 K=1,CP1 
ADFACI,K)=0,0 ^ 

DO 50 L=CP1,TN0C 

ALFACI,K)=Al!.FACl,R3+C0FrcI,t.)»BErA2CL,K3 
GO TO 80 

DO 60 I=CPl,T«OC 
DO 60 K=1,CP1 
ALFAtl,K3=0,0 
DO 90 K=1,CP1 
A6FACCT,KJs0,0 
DO 90 Li=CPl,TWOC 
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A|.FA(CT,K)sAl,FACCt,K3+COFF(CT,l,)»BEfA2(Ci,K3 

tiferUKrt ' ' * 

£N0 

SUBROUTiNE BCMULCJ) 

INTEGER C#CFl,TWttC,CT 
COMMON N,C,Cpl,T«OC,Cr 

^ r » ) /XCC/CC ( 9 , 9 ) /XR/P ( 8S , 9,5) 
RRRSCC ( i , C+1 ) ... 

CTCTsCCCCT.Cf) 

GO TO 100 

IFCA»5a,U-RRR3,LE.l,E-U6)GO TO 110 

DO '10 I=l f CT 

PC0fI,CPl5=BCI,CTlVCTCr 

DO 10 K=1,C 

PCJ,I,K)=RRR»BCI,K) 

DO 10 D=CF1,CT' 

10 P(J^lj^K)=F(J,I,K)-l«Ba,Dl*CCCL,K+C) 

HO DO 20 1=1, CT 

PCa/irCPl)s»Cl#CTJ»CTCT 
DO 20 K=l»C 
PCJ,I,K3s&(I,K) 

DO 20 D=CPt,CT 

20 PtUjj^I^KJsPCU^IjKJ+Ba^Dl^CCCL^K+C) 

100 IFCAB5(1.0-RRR),DE,1,E-06) 60 TO 120 

DO '30 1=1. CT 
PCJ£l,CPl5=BCl,Cf3*ClCT 
DO 30 K=l/C 

30 PCa^I|jK3sKRR#BCI,K) 

120 DO 40 1=1. CT 

PCJ,I,CT3=StI,CT)*CTCT 
DO 40 K=l,C 
40 Pta,l,KJ=8C3.,K3 

HETOKN 
END 

SUBROUTINE BAMUDCOZ. B.COFF) 

INTEGER C, CPI# TWdc.CT 
COMMON N,C,CP1 .TiiiUC.CT 
COMMON /XP/P(85,9,5J 

DIMENSION COFF(9,y).BC9,9) 

RRR=COFr(l,l) ‘ 

CTCTsCOrrtCT.CT) 

IFCAB5C1,0-RRR).OE.1.0E’"O6) GO TO 100 
DO 10 1=1. CT 

P(aZ,I^CP13=8(I.CT3»CTCT 

P?a2?I,R}=iRR»B(I.Kl+BCl,Cr)*COrF(CT,K3 
RETURN 

00 20 1=1. CT 

P(aZ.I.CPl)=BCI.CT3»CTCT 
00 '20 K«1.C 

|(OZ^I,N3=8tI.KJ4-BCI,CTl*COFrCCT.K3 
END 

SUBROUTINE MBAPDCUZl 
INTEGER C.CPlf TWOt.Cr 

COMMON /XP/BC§.93-^^ADrA/ALFAC9,S3/XP/PC85,9 
■DO 10 1=1, CT- • ' ■ ' 

PCjl?l£R3=*BCI.CT3!«=AOrAtCT,K3 

00 Iw W^lfC. \ 

10 PCJZ.I,K3aPtJZ.l.*^)-Ba.D3»ADFAlb.X) 

RETURN 
END 


10 

100 

20 


.53 
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SUBRUUTiriK MBApVfUZJ 
iNTt:Gt:K CfCFlyTWOC-CT 
COKHUN N,C,CP1,1’«UC|C'A’ 

uO 1 y I <• i ^ t. A ‘ _ . , . 

DO It# K=l»Ct'l 

PCaZ,I,K3=-BCl,CT)*At<rAtCT,K) 

DO 10 tiSlrC 

10 |t^2^1|‘^)=PCJZ#I,K)-Ba,03*A0fACIi,K) 

END 

SUBROOTINE AOMUDCUMl) 

INTEGER C,CPifTWOO,CT 

COMMON n,c,cp! ,T i«o£,eT 

COMMON /XA/AC9,»)/XU/a(47,9J/XVECl/VECl(9) 

RRR*ACljlA 

\fECHCT5sACCX,CTA*a(aMl»CT) 

SUMsO.O ' 

IFCAB5Cl.q-KRR),LE,l.fi-W6) GO TO 100 

DO 10 i»i,c 

S^EC1CI)-RHR»Q(JM1,I) 

VEC1CI+C)=0,0 

10 SUM=SUM +AICT,I)»0C0M1#A) 

VEClCCTlsVECKCTl+SOM 

RETURN 

100 DO 20 1=1, C 

VECl(I)=Q(aMl,I) 

VEC1CI+C)=0.0 

20 S0M=80M+ACCT,I)»QCaMl,IA 

VECl(CT)=V£ClCCT)'fSUM 
return 
END 

subroutine CQMUD(J) 

INTEGER C, CPI, TNOO.CT 
COMMON N-C,CPl,T»iaC,CT 

COMMON /XC6rF/CurE£9,9) /XBfiTAl/BETAlC93/XVECl/VEClC9) 
RRR rCOFFtl.CPi) ' " 

VECl (Cl }seorF(CT,CT) fBEf A1 (CT) 
IFCABSCl,a-ETA(a5nOE,l,E-06) GO TO 100 
IFCABS£1.0-KRR3*DE.1,E*06) go to 50 
DO -10 1=1, C 

VECitI)=RKR»BETAUI + C) 

VEC1CI+C)=0,0 

VECl CCT3=VEC1CCT)+C0FF£CT,I+C)*»ETA1(I+C) 

DO 10 L.=CP1,TW0C ■ - 

10 VECICI + C)=VECI (l4-C)+C0FF£ltC,Ii)4BErAlCl.) 

RETURN ' ■ • 

50 DO 20 1 = 1, C 

VECI(I)=BETA1(I+C) 

VECI(I+C)=0,0 

FECI tCT)=VECl ( Cf ) +COFF£CT, I+C3 »BETA1 C I+C) 

DO 20 D=CPl,TWOC 

20 FECHI+C)=VEC1(I4-C)+CDFF(I+C,L)»BEXA1CD) 

RETURN ^ ' ' ' 

100 IFCABS£l,0«-KRR3,DE,l,E-05) GO TO 150 

DO *30 1=4, C 

FECI tI>=R^R»BETAl tl+CJ 
FECI (14-0=0*0 

30 FEClCCT)=VECl(CT3+C0FrCCT,I+C)*BETAl(I+O 

RETURN - 

150 DO 40 1=1, C ,, 

FECl(I)=iETAiCl4-C) 

FECI (I tC) =0.0 

40 VEClCCT)=FEClCCT)+C0FrCCT,I+C)4BEfAl(I+C) 

RETURN ' - 

END 

SUBROUTINE BFMUDCa) 

INTEGER C,CPl,TWOC,CT 
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COMMON N,C,Cpl,!rWUC/Cr 

I/O Id 

acj#ij=o,u 

DO 10 K=lrCl' 

END 

subroutine OMPVXCU) 

INTEGER CfCPl.TNOC'cT 

COMMON N-C,CpI,'r«OC,CT 

COMMON /XP/PC85,9,5)/XQ/a(47,9) 

DO 10 1=1 fCT 
DO 10 K=lfCFl 

10 0(J^^I)=UCJ,X)-P(J,I,K)fU(J+l,!(*C) 

END 

SUBROUTINE aMPiiXca,ap,au) 

INTEGER C,CPl»TWOC^Cf 

COMMON N-C,CPl,TWCJC»Cf 

COMMON /XP/PC85,9,5}/XQ/Q(47,9) 

DO Id ' 

QCJ|I)=ac5^1)-PCaP,l,CPl)»QCD«rCf5 
10 Q(a,i}p«c5,i)-p{aP,i,K>»Q(jQ,K) 

RETURN - • ' 

END ■ 

SUBROUTINE INVPRflBJ 
INTEGER C,CT 
COMMON C,CI,N 

DIMENSION 8CCT,CT),BSt4#4),8SlC4,4) ,PA(4),PBC4) 

DO 1 1 = 1, C ' '^ * ■ • - . . 

PA(i)=u/eti,i) 

CONTINUE 
DO 2 1=1, C 
DO 2 J=l,C 

BStI, J)=PA(1)»BCI#J+C) 

J CONTINUE 

DO 3 I=1,C 
DO 3 Jsl.C 
DO 3 ia=i,c 

BSlCI,J3=B51(I,J)+BlH*C,ia)»BS(ia,a) 

i CONTINUE ^ 

DO 4 1=1, C 
00 4 J=X,C 

BSlCI,J)=8(I+C,a+C)-BSltI,vJ) 

I CONTINUE 

CADE MATINC851,C, XXX, 0, DETERM) 

DO b 1=1, C- 
DO 5 Jsl.C 
B(I'»'C,J+C)=BSl(I,D) 
asici,a)=Bci+c,-j)»pA(j) 

eti,u+c)=o, 

BCi4-c,a)=o, 

5 CONTINUE 

DO 6 1=1, C 
00 & J=l,C 

B?i£j+c)=6ci,a+c)-Bsu,iJ?»Bcia+c> j+c) 
b continue ' 

DO 7 I=1,C 
DO 7 u=i;c 
8S(I,J)a6, 

BClJc^J)=Aci+C,J)-BU+C,C+IU3«BSlCIJ,a) 
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DO "9 1=1, c 
BSCI.I)=^A(i) 

pAci5=fi. 

DO '9 J=J, ,C 

DO ID 1=1, C 
DO 10 11=1, C 

PMI)=PAa5+BCI,II)»aCll,Ct)fBCl,II+C)tB(lI+C,CT) 

wuntiwuc* - - 

00 11 1=1, C 

^ggjB8j|CCl,l)»PAU)+B(CT,l+C)»=PBCl) 

DO 1 2 I”! , C 

8Cl,CT)=-PA(l)»BC0f,CI} 

8 ( j +c , cf ) =-PB ( I ) a c cr , ci‘ ) 

PACX)=0. 

PBCI)=0, 

CONflNUE 
DO 13 J=1,C 
DO 13 11=1 C 

PA(J3=PA(jUBtCT,II)»BClI,J)4-B(CT,Il+C)»BC2I+C,a) 

PBCJJ=PBCDHBCCT;iIJ»BCli;a+C)+8(^T,II+C)»BCII+C,i 

B?CX?I)=-S?CT,CT)»PA(1) 

DO 15 1=1, C 
DO 15 a=l,C 

B(I,D3=BCI,a)-BCI,Cf )»PACJ) 
BCI,J+C3=BU,J+C)«Ba,C'l*)*PB(a3 
8CI+C,J}=BCI+C,a)-BCl+C,CT)*PAca3 
8CI+C,a+C3=BCl+C,vJ+C3-BtI+C,CT3»PB(a) 

8S1CI,J3=0, 

COMTlflUE 

BB=0. 

DO 22 1=1, C 
PBC1J=0, 

COKTINUK 

RETURN 

END 


a+c) 


SUBROUTINE MATIN C A, N,B,«, DETERM) 

A=C0-EFriCieNT OF OKOEr N 
B=VECTOK or ORDER N 

M=IF M IS set to 5iEK0,0N0y INVERSEIS COMPUTED 
DETERMsVAOUE of DETERMENENT RETURNED 

DIMENSION At N,N3,BCN,1) »IPIV0T( TO 3, INDEXC70,2),0Tl703 

EOUIVAoENCE CIROW,0k6n) , CICODUM, JC0DUM3 , CAMAX,T,SWAP3 

INITIAOIEATAON 

DETER«=l,0 

DO 20 0=1, N 

IPIVOTCU3=0 

SEARCH FOR PIVOT ELEMENT 
DO 550 1=1, N 
AHAX=U.0 
DO 105 U=1,!N 

IFClPIVOTtai-l 360, 105,50 
DO 100 K=1,N 

lFtlPIVOTtK3-l)«0, 100,740 
IFCAMAX-ABSIA(vI,k 53 J85,100,100 
IROWsJ 
ICOLUM=K 
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100 

IU5 

C 

140 

200 

210 


250 

260 


220 

C 

360 

370 

380 

400 


450 

460 

500 

550 

C. 


630 


705 

710 


11 

12 

991 

740 

C 


ftMAX=AB6CA(U,K)) 

CONTINUE ■ ' 

CONTINUE 

IPIVOTtICOLOM)=IPlVOI(lC0liUMJ+l 

INTRKCHANOE rows fO POT PlyOT VECTORS ONOIAGONAL 
IFClR0w-lC00UMJl40-2e0;i40 uwyxauuflAi, 

DETERMr-DElERM ' - ' 

DO 200 0=1 -K 
SWAPxAClKO^^ 

■ ■ ■ ■ IrOJ 


SWAP=ACIKOW,OJ 

ACIHUW,L,)=AtICCILU«, 

AClCOOUM,li)aSWAP 

T » # lyi % n I: rut 


IFC«}260;260,210 
DO- 250 L=i,R 
SWAPsSClROWfU) 

8ClROW--li)=B(IC01iU«,li5 
B(ICOliUM,0)*SWAP 
INDEXCI,iJ=lROW 
XN0EXC I »2) sICOLUM 
DIVIDE ^IVOT ROW PIVOR EDEMg«T 
PIVOX=A(ILmU«,ICOLUM) 

DTCI)»PIVUT 
A(:iCOLUM,iCO0UM?al,O 
DO 220 lisl.N 

A(ICO0U«,ti5=ACIC0LUM,0) /PIVOT 
REDUCE NOW PIVOT KO»S 
IFCA)380,380,360 
DO-370 L«1,H 

B(ICODUM,L)=a(ICOOUM,l.) /PIVOT 
DO 550 1*1 = 1, M 
IFCI*l-ICOIi0M)4OO,55O,4OO 
TsACDltlCOLUfi) 

ACIil,ICODOM>=0,0 
DO 450 0=1, N 

AC01,O)=At0l,L3-AlICOOUN,O)»r 
IFCm5550,550,460 
DO -500 0=1, m' 

BCOl,L)=BlLl,03«BUCOOUM,0)*f 
CONTINUE 

INTERCHANGE THE COOUMNS 
DO 710 1=1, N 
0=N+l»l 

IFCINDEXC0,13-INDEXt 0,23 3630,710,630 
JROW=INOEXtO,l3 ' 

3COOUM=INUEXtO,2) 

DO 705 K=l,N 
SWAP=A(K,JftOW) 

A(K,JRQW3=A1K, JCOOUM) 

ACK,aC0LiUM3=SWAP 
CONTINUE 
CONTINUE 
DO 11 K=1,N 

iF(IPIVOTtK3,NE,l3GO TO 12 
CONTINUE ' 

RETURN 

WRIfE(22,99a3 

F0RMATC/3OX, ^MATRIX IS SiNGUOAR* 

RETURN 
END 

#»:(! f! 3|£ If nt K! ♦ W f »***» f » f »*»♦♦ m **»»#»»»*♦*»♦♦» )K ♦» w f w f * f »» »t i 
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